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Abstract 

Recent experimental data have estabhshed some of the basic features of mul- 
tiple interactions in hadron-hadron collisions. The emphasis is therefore now 
shifting, to one of exploring more detailed aspects. Starting from a brief re- 
view of the current situation, a next-generation model is developed, wherein a 
detailed account is given of correlated flavour, colour, longitudinal and trans- 
verse momentum distributions, encompassing both the partons initiating per- 
turbative interactions and the partons left in the beam remnants. Some of 
the main features are illustrated for the Tevatron and the LHC. 
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Figure 1: Schematic illustration of an event with two 2 — > 2 perturbative interactions. 

1 Introduction 

The physics of high-energy hadron-hadron interactions has become a topic of increasing 
interest in recent years. With the Tevatron Run II well under way and with the startup of 
the LHC drawing closer, huge data samples are becoming available that will challenge our 
current understanding of this physics. From the point of view of QCD, many interesting 
questions remain to be answered, and we shall take up some of these in detail below. 
Moreover, for new physics searches and precision measurements, it is important that these 
questions can be given meaningful and trustworthy answers, since ever-present yet poorly- 
understood aspects of QCD can have a significant impact. 

Much of the complexity involved in describing these phenomena — specifically the 
underlying event and minimum-bias collisions — derives from the composite nature of 
hadrons; we are dealing with objects which possess a rich internal structure that is not 
calculable from perturbation theory. This, however, does not imply that the physics of the 
underlying event as such has to be an inherently non-perturbative quagmire. 

Viewing hadrons as 'bunches' of incoming partons, it is apparent that when two hadrons 
collide it is possible that several distinct pairs of partons collide with each other, as depicted 
in Fig. 1. Thus multiple interactions (also known as multiple scatterings) in hadronic col- 
lisions is a phenomenon which is a direct consequence of the composite nature of hadrons 
and which must exist, at some level. In fact, by extending simple perturbation theory 
to rather low pj_ values, though still some distance above Aqcd? most inelastic events in 
high-energy hadronic collisions are guaranteed to contain several perturbatively calculable 
interactions [1]. Furthermore, such interactions — even when soft — can be highly impor- 
tant, causing non-trivial changes to the colour topology of the colliding system as a whole, 
with potentially drastic consequences for the particle multiplicity in the final state. 

Nevertheless, traditionally the exploration of multiple interactions has not attracted 
much interest. For studies concentrating on high-p^ jets, perturbative QCD emission is 
a more important source of multijets than separate multiple interactions. The underlying 
event, on the other hand, has in this context often been viewed as a mess of soft QCD in- 
teractions, that cannot be described from first principles but is better simply parametrized. 

However, such parametrizations, even while reasonably successful in describing the av- 
erage underlying activity, are not sophisticated enough to adequately describe correlations 
and fluctuations. This relates for instance to jet profiles and jet pedestals, and to systematic 
as well as random shifts in jet energies. The lack of sophistication implies that, even when 
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tuned to describe a few distributions well, they could not be trusted for extrapolations 
beyond the fit region. Hence, a sound understanding of multiple interactions is prerequisite 
for precision physics involving jets and/or the underlying event. 

It is interesting to note that this can also impact physics studies in areas well beyond 
the conventional QCD ones. As an example, consider the search for a Higgs particle in the 
h° 77 channel at the LHC, where the Higgs mass resolution at high luminosity depends 
on picking the correct vertex between maybe 30 different pp events. If the event that 
contained the Higgs is special, by typically containing more charged particles (in general 
or in some regions of phase space), we would want to use that information to clean up the 
signal [2]. 

The crucial leap of imagination is to postulate that all particle production in inelastic 
hadronic collisions derives from the multiple-interactions mechanism. This is not to say 
that many nonperturbative and poorly known phenomena will not force their entrance on 
the stage, in going from the perturbative interactions to the observable hadrons, but the 
starting point is perturbative. If correct, this hypothesis implies that the typical Tevatron 
hadron-hadron collision may contain something like 2-6 interactions, and the LHC one 
4-10. 

A few models based on this picture were presented several years ago [1] , and compared 
with the data then available. Though these models may still be tuned to give a reasonable 
description of the underlying event at various collider energies, several shortcuts had to be 
taken, particularly in the description of the nonperturbative aspects alluded to above. For 
instance, it was not possible to consider beam remnants with more than one valence quark 
kicked out. 

The increased interest and the new data now prompts us to develop a more realistic 
framework for multiple interactions than the one in ref. [1], while making use of many of 
the same underlying ideas. This may not necessarily result in significant improvements for 
fits made to only a few distributions at a time, but we hope it will enhance our ability to 
simultaneously describe many different observations inside one framework, thereby improv- 
ing the confidence with which we can make extrapolations from known measurements to 
new distributions and to higher energies. 

One of the building blocks for the new model comes from our recent study of baryon- 
number- violating processes [3] . We then had to address the hadronization of colour topolo- 
gies of the same kind as found in baryons. Specifically, as an extension of the standard 
Lund string fragmentation framework [4], we explored the concept of a junction in which 
three string pieces meet, with a quark at the end of each string. This also opens the way 
to a more realistic description of multiple interaction events. 

The resulting improvements, relative to the framework in [1], are especially notable 
in the description of the structure of the incoming hadrons, i.e. how flavours, colours, 
transverse and longitudinal momenta are correlated between all the partons involved, both 
those that undergo interactions and those that are left behind in the remnants. (Brief 
descriptions of some of these aspects can also be found in [5].) To give one specific example, 
we introduce parton densities that are modified according to the flavours already affected 
by interactions. 

Clearly, the model we present here is not the final word. For instance, we defer for 
the future any discussions of processes involving photons in the initial state and whether 
and how diffractive topologies could arise naturally from several interactions with a net 
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colour singlet exchange. More generally, the whole issue of colour correlations will require 
further studies. The model also allows some options in a few places. A reasonable range 
of possibilities can then be explored, and (eventually) experimental tests should teach us 
more about the course taken by nature. 

This article is organized as follows. In Section 2 we give an introduction to multiple 
interactions in general, to the existing multiple interactions machinery, to other theoretical 
models, and to experimental data of relevance. Sections 3-5 then describe the improvements 
introduced by the current study: 3 a new option for impact-parameter dependence, 4 the 
main work on flavour and momentum space correlations, and 5 the very difficult topics 
of colour correlations and junction hadronization. Finally Section 6 provides some further 
examples of the resulting physics distributions and important tests, while Section 7 contains 
a summary and outlook. 

2 Multiple Interactions Minireview 
2.1 Basic Concepts 

The cross section for QCD hard 2^2 processes, as a function of the scale, is given by 

^ = J]|dxi|dx2|dtV,(xi,g^)/,(x2,g^)^^5(pi-y) , (1) 

where s = X1X2S. The jet cross section is twice as large, crjct = 2crint, since each interaction 
gives rise to two jets, to first approximation. In the following, we will always refer to the 
interaction rather than the jet cross section, unless otherwise specified. We will also assume 
that the 'hardness' of processes is given by the p± scale of the interaction, i.e. = p\. 

The cross section for QCD 2 — 2 processes, the sum of qq' — > qq', qq q'q', qq — >• gg, 
qg — s> qg, gg — > gg and gg qq, is dominated by t-channel gluon exchange contributions. 
(This justifies the use of 'interaction' and 'scattering' as almost synonymous.) In the \t\ <^ s 
limit, where p\ = tu/s ~ \t\, quark and gluon interactions just differ by the colour factors, 
so approximately we may write 



where 



da _ Snalipl) 



(3) 



and 



dp]_ 9p\ 

F{x, g2) = 5^ (x q{x, g2) + xqix, Q^)) + lxg{x, Q^) . (4) 

Q 

Thus, for constant and neglecting the x integrals, the integrated cross section above 
some p_Lmin is divergent in the limit p± ^ 0: 

da 1 

0-int(p±min) = / dp_L OC — . (5) 
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Figure 2: The integrated interaction cross section crint above p_Lmin for the Tevatron, 
with 1.8 TeV pp colhsions, and the LHC, with 14 TeV pp ones. For comparison, the 
flat hnes represent the respective total cross section. 



A numerical illustration of this divergence is given in Fig. 2. Note that the actual fall- 
off is everywhere steeper than l/p\- We have here used full 2^2 matrix elements and 
the CTEQ 5L parton density parametrizations [6], which are valid for Q > 1.1 GeV and 
X > 10~^; therefore results at the lowest pj_ values are not to be taken too literally. For the 
studies in this article we are basing ourselves on leading-order cross sections and parton 
densities, with nontrivial higher-order corrections only approximately taken into account by 
the addition of parton showers. Nevertheless, the trend is quite clear, with an integrated 
cross section that exceeds the total pp/pp cross section atot (in the parametrization of 
ref. [7]) for p^min of the order of a few GeV. As already mentioned in the introduction, this 
is well above Aqcd? so one cannot postulate a breakdown of perturbation theory in the 
conventional sense. 

The resolution of the (Tint > o'tot paradox probably comes in two steps. 

Firstly, the interaction cross section is an inclusive number. Thus, if an event contains 
two interactions it counts twice in o"int but only once in atot, and so on for higher multi- 
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plicities. Thereby we may identify (n)(p_Lmin) = crint(p±min)/crtot with the average number 
of interactions above p_Lmin per event, and that number may well be above unity. 

One of the problems we will consider further in this article is that this simple cal- 
culation of {n){p±min) does not take into account energy-momentum conservation effects. 
Specifically, the average s of a scattering decreases slower with p_Lmin than the number of 
interactions increases, so naively the total amount of scattered partonic energy becomes 
infinite. Thus corrections reduce the (n)(p_Lmin) number, but not sufficiently strongly: one 
is lead to a picture with too little of the incoming energy remaining in the small-angle beam 
jet region [1]. 

Secondly, a more credible reason for taming the rise of (?^)(p_Lmin) is that the incoming 
hadrons are colour singlet objects. Therefore, when the p± of an exchanged gluon is made 
small and the transverse wavelength correspondingly large, the gluon can no longer resolve 
the individual colour charges, and the effective coupling is decreased. Note that perturbative 
QCD calculations are always performed assuming free incoming and outgoing quark and 
gluon states, rather than partons inside hadrons, and thus do not address this kind of 
nonperturbative screening effects. 

A naive estimate of an effective lower cutoff would be 

h 0.2 GeV-fm ^ , ^ , 

P±min ~ — ^ -— ^ 0.3 GeV ~ Aqcd , 6 

Tp 0.7 fm 

but this again appears too low. The proton radius Vp has to be replaced by the typical colour 
screening distance d, i.e. the average size of the region within which the net compensation of 
a given colour charge occurs. This number is not known from first principles, so effectively 
one is forced to introduce some kind of cutoff parameter, which can then just as well be 
put in transverse momentum space. The simplest choice is to introduce a step function 
9(j)± — p±min), such that the perturbative cross section completely vanishes below the p_Lmm 
scale. A more realistic alternative is to note that the jet cross section is divergent like 
'^s(pi)/Pj_5 that therefore a factor 

would smoothly regularize the divergences, now with p±Q as the free parameter to be tuned 
to data. Empirically the two procedures give similar numbers, p_Lmin ~ P±o, and both of 
the order of 2 GeV. 

The parameters do not have to be energy- independent, however. Higher energies imply 
that parton densities can be probed at smaller x values, where the number of partons 
rapidly increases. Partons then become closer packed and the colour screening distance d 
decreases. Just like the small-x rise goes like some power of x one could therefore expect 
the energy dependence of p_Lmin and p±o to go like some power of CM energy. Explicit toy 
simulations [8] lend some credence to such an ansatz, although with large uncertainties. 
Alternatively, one could let the cutoff increase with decreasing x; this would lead to a 
similar phenomenology since larger energies probe smaller x values. 

2.2 Our Existing Models 

The models developed in ref. [1] have been implemented and are available in the Pythia 
event generator. They form the starting point for the refinements we will discuss further 
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on, so we here review some of the main features. 

The approach is not intended to cover elastic or diffractive physics, so the crmt{p±ram, s) or 
c"int(p±0) s) interaction cross section is distributed among the crnd(s) nondiffractive inelastic 
one [7, 10]. The average number of interactions per such event is then the ratio (n) = 
o"int/o"nd- As a starting point we will assume that all hadron collisions are equivalent, i.e. 
without an impact parameter dependence, and that the different parton-parton interactions 
take place independently of each other. The number of interactions per event is then 
distributed according to a Poisson distribution with mean (n), Vn = (n)'"exp(— 

One (not used) approach would be, for each new event, to pick the actual number 
of interactions n according to the Poissonian, and select the n pj_ values independently 
according to eq. (1). One disadvantage is that this does not take into account correlations, 
even such basic ones as energy-momentum conservation: the sum of interaction energies 
may well exceed the total CM energy. 

In an event with several interactions, it is therefore convenient to impose an ordering. 
The logical choice is to arrange the scatterings in falling sequence of p± values. The 'first' 
scattering is thus the hardest one, with the 'subsequent' ('second', 'third', etc.) successively 
softer. This terminology is not primarily related to any picture in physical time although, 
by the uncertainty relation, large momentum transfers implies short timescales. When 
averaging over all configurations of soft partons, one should effectively obtain the standard 
QCD phenomenology for a hard scattering, e.g. in terms of parton distributions. Correlation 
effects, known or estimated, can be introduced in the choice of subsequent scatterings, given 
that the 'preceding' (harder) ones are already known. This will be developed further in 
Section 4. 

The generation of a sequence y/s/2 > p±i > p±2 > ■ ■ ■ > p±n > P±min now becomes one 
of determining p±i from a known p±i-i, according to the probability distribution 



dV 1 da 
dp±i 0-nd dp± 



P^'-' 1 d(T , , 

CTnd dp. 



(8) 



The exponential expression is the 'form factor' from the requirement that no interactions 
occur between p±i-i and p±i, cf. radioactive decays or the Sudakov form factor [11] of 
parton showers. 

When used with the standard differential cross section dcr/dp_L, eq. (8) gives the same 
Poisson distribution as above. This time n is not known beforehand, but is defined by the 
termination of the iterative procedure. Now, however, da/dp± can be modified to take into 
account the effects of the i — 1 preceding interactions. Specifically, parton distributions are 
not evaluated at Xi for the i'th scattered parton from a hadron, but at the rescaled value 

= . X. . (9) 



Xj 



so that it becomes impossible to scatter more energy than initially available in the incoming 
beam. This also dynamically suppresses the high-multiplicity tail of the Poissonian and 
thereby reduces the average number of interactions. 

In a fraction of the events studied, there will be no hard scattering at all above p_Lmin- 
Such events are associated with nonperturbative low-p_|_ physics, and are simulated by ex- 
changing a very soft gluon between the two colliding hadrons, making the hadron remnants 
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colour-octet objects rather than colour- singlet ones. If only valence quarks are considered, 
the colour-octet state of a baryon can be decomposed into a colour triplet quark and an 
antitriplet diquark. In a baryon-baryon collision, one would then obtain a two-string pic- 
ture, with each string stretched from the quark of one baryon to the diquark of the other. 
A baryon-antibaryon collision would give one string between a quark and an antiquark and 
another one between a diquark and an antidiquark. 

In a hard interaction, the number of possible string drawings are many more, and 
the overall situation can become quite complex. In the studies preceding this work, several 
simplifications were made. The hardest interaction was selected with full freedom of flavour 
choice and colour topology, but for the subsequent ones only three simple recipes were 
available: 

• Interactions of the gg gg type, with the two gluons in a colour-singlet state, such 
that a double string is stretched directly between the two outgoing gluons, decoupled 
from the rest of the system. 

• Interactions gg — gg, but colour correlations assumed to be such that each of the 
gluons is connected to one of the strings 'already' present. Among the different 
possibilities of connecting the colours of the gluons, the one which minimizes the total 
increase in string length is chosen. This is in contrast to the previous alternative, 
which roughly corresponds to a maximization (within reason) of the extra string 
length. 

• Interactions gg — > qq, with the final pair again in a colour- singlet state, such that a 
single string is stretched between the outgoing q and q. 

The three possibilities can be combined in suitable fractions. 

Many further approximations were also required in the old framework, e.g. the addition 
of initial- and final-state parton showers was feasible only for the hardest interaction, and 
we address several of those in the following. 

The model also includes several options for the impact-parameter dependence. This 
offers an additional element of variability: central collisions on the average will contain 
more interactions than peripheral ones. Even if a Poisson distribution in the number of 
interactions would be assumed for each impact parameter separately, the net result would 
be a broader-than-Poisson distribution. The amount of broadening can be 'tuned' by the 
choice of impact-parameter profile. We discuss this further in section 3, where a new set of 
profiles is studied. 

The above framework was originally formulated for pp/pp collisions, but has also been 
extended to 7p and 77 interactions [12]. In these latter processes, however, the nature of 
the photon needs to be modelled in detail, and this introduces many further uncertainties. 
A study of such aspects is beyond the scope of the current article. 

2.3 Other Models 

While the models of ref. [1] may well be the ones most frequently studied, owing to their 
implementation in Pythia [9], a number of other models also exist. Many of the basic 
concepts have also been studied separately. We here give a few examples, without any 
claim of completeness. 

In Dual Topological Unitarization (DTU) language [13], and the Dual Parton Model 
based on it [14], or other similar techniques [15], inelastic events are understood in terms 
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of cut pomerons [16]. Translated into modern terminology, each cut pomeron corresponds 
to the exchange of a soft gluon, which results in two 'strings' being drawn between the two 
beam remnants. Uncut pomerons give virtual corrections that help preserve unitarity. A 
variable number of cut pomerons are allowed. This approach has been the basis for the 
simulation of underlying events in ISAJET [17], and was the starting point for Dtujet [18]. 
However, note that cut pomerons were originally viewed as purely soft objects, and so 
did not generate any transverse momentum, unlike the multiple interactions considered in 
this article. In Dtujet and its Phojet [19] and Dpmjet [20] relatives, however, also 
hard interactions have been included, so that the picture now is one of both hard and 
soft pomerons, ideally with a smooth transition between the two. Since the three related 
programs make use of the Pythia hadronization description, the differences relative to our 
scenarios is more a matter of details (but "the devil is in the details") than of any basic 
incompatibility. 

The possibility of observing two separate hard interactions has been proposed since long 
[21], and from that has also developed a line of studies on the physics framework for having 
several hard interactions [22], also involving e.g. electroweak processes [23]. Again this is 
similar to what we do, except that lower p± values and the transition to nonperturbative 
physics are not normally emphasized. 

The possibility of multiple interactions has also been implicit [24] or explicit [25] in 
many calculations of total cross sections for hadron-hadron, hadron-7 and 77 events. The 
increase of a^^t with CM energy is here directly driving an increase also of (Ttot; that the 
latter is rising slower than the former comes out of an eikonalization procedure that implies 
also an increasing (riint). 

Multiple interactions require an ansatz for the structure of the incoming beams, i.e. 
correlations between the constituent partons. Some of these issues have been studied, e.g. 
with respect to longitudinal momenta [14,26], colours [27] or impact parameter [28], but 
very little of this has been tested experimentally. Dense-packing of partons could become 
an issue [15], of unknown importance, but not necessarily a major one [29]. 

The Herwig [30] event generator does not contain any physics simulation of multiple 
interactions. Instead a parametrization procedure originally suggested by UA5 [31] is used, 
without any underlying physics scenario. It thus does parametrize multiplicity and rapidity 
distributions, but does not contain any minijet activity in the underlying event. The add- 
on Jimmy package [32] offers a multiple-interaction component (both for pp, 7p and 77 
events), which has recently been extended to include also a model of soft interactions [33]. 

The introduction of unintegrated parton densities, as used in the BFKL/CCFM/LDC 
approaches to initial-state radiation [34-36], allows the possibility to replace our p±min/p±o 
cutoff by parton densities that explicitly vanish in the p± ^ limit [37]. This opens up 
the possibility of an alternative implementation of multiple interactions [38]. 

In heavy-ion collisions the multiple interactions rate can become huge [39]. For small 
impact parameters, a major fraction of the energy of the two incoming nuclei is carried 
by partons undergoing perturbative interactions, and which therefore define a 'resolved' 
partonic content. This suggests a mechanism for the construction of an 'initial state' for 
the continued formation and thermalization (or not) of a quark-gluon plasma. 
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2.4 Experimental Tests 

Experimental input to the understanding of multiple interactions comes in essentially three 
categories: direct observation of double parton scattering, event properties that directly 
and strongly correlate with multiple interactions, and event properties that do not point to 
multiple interactions by themselves but still constrain multiple interactions models. 

If an event contains two uncorrelated 2 — *■ 2 interactions, we expect to find four jets, 
grouped in two pairs that each internally have roughly opposite and compensating trans- 
verse momenta, and where the relative azimuthal angle between the scattering planes is 
isotropic. Neither of these properties are expected in a 2 ^ 4 event, where two of the 
partons can be thought of as bremsstrahlung off a basic 2^2 process. The problem is 
that 2—^4 processes win out at large p±, so there is a delicate balance between having 
large enough jet p± that the jets can be well measured and still not so large that the signal 
drowns. 

When the p±min of the jets is sufficiently large that exp(— (n)) ~ 1, Poisson statistics 
implies that V2 = 'Pf/2, where Vi is the probability to have i interactions. Traditionally 
this is rewritten as 



where the ratio (Xnd/c'cfr gauges deviations from the Poissonian ansatz. Values above unity, 
i.e. (Jeff < cTjid, arise naturally in models with variable impact parameter. 

The first observation of double parton scattering is by AFS [40]. The subsequent UA2 
study [41] ends up quoting an upper limit, but has a best fit that requires them. A CDF 
study [42] also found them. These experiments all had to contend with limited statistics 
and uncertain background estimates. The strongest signal has been obtained with a CDF 
study involving three jets and a hard photon [43]; here a2 = (TAcrB/crcs, without a factor 
1/2, since the two 2^2 processes A and B are inequivalent. In all cases, including the 
UA2 best fit, cXeflf comes out smaller than Und; typically the double parton scattering cross 
section is a factor of three to four larger than the Poissonian prediction. For instance, the 
CDF number is o"efr = 14.5 ± 1.7l2'3 More recently, ZEUS has observed a signal in 7p 
events [44]. The DO four-jet study shows the need to include multiple interactions, but this 
is not quantified [45]. 

So far, no direct tests of triple or more parton scattering exist. However, the UAl 
minijet rates [46], going up to 5 jets, are difficult to understand without such events. 

Tests involving jets at reasonably large p± values do not probe the total rate of multiple 
interactions, since the bulk of the interactions occur at so small p± values that they cannot 
be identified as separate jets. By the way colours are drawn across the event, soft partons 
can drive the production of particles quite out of proportion to the p_i_ values involved, 
however. The multiplicity distribution of multiple interactions thereby strongly influences 
the multiphcity distribution of charged hadrons, rich- 

A notable aspect here is that the measured rich distribution, when expressed in the KNO 
variable z = rich/inch) [47], is getting broader with increasing CM energy [48,49]. This is 
contrary to the essentially Poissonian hadronization mechanism of the string model, where 
the KNO distribution becomes narrower. As an example, consider the UA5 measurements 
at 900 GeV [48], where {rich) = 35.6 and a{nch) = 19.5, while the Poissonian prediction 
would be cr(nch) = ^/(nd^ = 6.0. It is possible to derive approximate KNO scaling in 
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e"''e~ annihilation [50] , but this then rests on having a perturbative shower that involves 
the whole CM energy. However, allowing for at most one interaction in pp events and 
assuming that hadronization is universal (so it can be tuned to e+e~ data), there is no 
(known) way to accommodate the experimental multiplicity distributions, neither the rapid 
increase with energy of the average nor the large width. Either hadronization is very 
different in hadronic events from e"'"e~ ones, or one must accept multiple interactions as 
a reality. (Unfortunately it is difficult to test the 'hadronization universality' hypothesis 
completely separated from the multiple interactions and other assumptions. To give two 
examples, the relative particle flavour composition appears to be almost but not quite 
universal [51], and low-mass diffractive events display 'string-like' flavour correlations [52].) 

Further support is provided by the study of forward-backward multiplicity correlations. 
For instance, UA5 and E735 define a 'forward' np and a 'backward' ub multiplicity in 
pseudorapidity windows of one unit each, separated by a Ar] variable-width gap in the 
middle [53]. A forward-backward correlation strength is now defined by 



where the last equality holds for a symmetric r] distribution, i.e. for pp/pp but not for 
7p. Measurements give a positive and surprisingly large b, also for At] of several units of 
rapidity. So it appears that there is some global quantity, different for each event, that 
strongly influences particle production in the full phase space. Again known fragmentation 
mechanisms are too local, and effects of a single hard interaction not strong enough. But 
the number of multiple interactions is indeed a global quantity of the desired kind, and 
multiple-interaction models can describe the data quite well. 

It is a matter of taste which evidence is valued highest. The direct observation of double 
parton scattering is easily recognized as evidence for the multiple-interactions concept, but 
only affects a tiny fraction of the cross section. By comparison, the broadening multiplicity 
distribution and the strong forward-backward correlations offer more indirect evidence, but 
ones strongly suggesting that the bulk of events have several interactions. We are not aware 
of any realistic alternative explanations for either of the observables. 

Another interesting phenomenon is the pedestal effect: events with high-p^ jets on the 
average contain more underlying activity than minimum-bias ones, also well away from jets 
themselves. It has been observed by several collaborations, like UAl [54], CDF [55,56] and 
HI [57]. When the jet energy is varied from next-to-nothing to the highest possible, the 
underlying activity initially increases, but then flattens out around p_Ljct = 10 GeV (details 
depending on the jet algorithm used and the CM energy). This fits very nicely with an 
impact-parameter-dependent multiple-interactions scenario: the presence of a higher p± 
scale biases events towards a smaller impact parameter and thereby a higher additional 
activity, but once crint(p±jct) ^ CTnd the bias saturates [1]. The height of the pedestal 
depends on the form of the overlap function 0{b), and can thus be adjusted, while the 
P±}ct at which saturation occurs is rather model-insensitive, and in good agreement with 
the data. 

The presence of pedestals also affects all measurements of jet profiles [55,59]. It can 
lead to seemingly broader jets, when the full underlying event cannot be subtracted, and 
enrich the jet substructure, when a multiple-interactions jet is mistaken for radiation off 
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the hard subprocess. It can also affect (anti) correlations inside a jet and with respect to 
the rest of the event [57]. 

Many further observables influence the modeling and understanding of multiple inter- 
actions, without having an immediate interpretation in those terms. A notable example 
here is the {p±){nch) distribution, i.e. how the average transverse momentum of charged 
particles varies as a function of the total charged multiplicity. The observed increasing 
trend [58] is consistent with multiple interactions: large multiplicity implies many interac- 
tions and therefore more perturbatively generated p± to be shared between the hadrons. 
For it to work, however, each new interaction should add proportionately less to the total 
rich than to the total p±. Whether this is the case strongly depends on the colour connec- 
tions between the interactions, i.e. whether strings tend to connect nearest neighbours in 
momentum space or run criss-cross in the event. A rising trend can easily be obtained, but 
it is a major challenge to get the quantitative behaviour right, as we shall see. 

Finally, one should mention that global fits to hadron collider data [60-64] clearly point 
to the importance of a correct understanding of multiple interactions, and constrains mod- 
els down to rather fine details. This brings together many of the aspects raised above, 
plus some further ones. A convenient reference for our continued discussion is Tune A, 
produced by R.D. Field, which is known to describe a large set of CDF minimum bias 
and jet data [61]. Relative to the defaults of the old scenario, it assumes p±o = 2.0 GeV 
(PARP(82)=2 . 0) at the reference energy 1.8 TeV (PARP(89)=1800 . 0), with an energy rescal- 
ing proportional to E^m (PARP(90)=0 . 25). It is based on a double Gaussian matter distri- 
bution (MSTP(82)=4), with half of the matter (PARP(83)=0 . 5) in a central core of radius 
40% of the rest (PARP(84)=0.4). Almost all of the subsequent interactions are assumed 
to be of the type gg — > gg with minimal string length (PARP(85)=0.9, PARP(86)=0 . 95). 
Finally the matching of the initial-state showers to the hard scattering is done at a scale 

Qshower = 4pihard (PARP (67) =4 . O) . 

The above parameter set is sensible, within the framework of the model, although by 
no means obvious. The matter distribution is intermediate between the extremes already 
considered in [1] , while the string drawing is more biased towards small string lengths than 
foreseen there. The p±Q energy dependence is steeper than previously used, but in a sensible 
range, as follows. In reggeon theory, a Pomeron intercept of 1+e implies a total cross section 
rising like s*^, and a small-x gluon density like xg{x) oc x"*^ (at small Q"^). A p_Lo rising (at 
most) like s'^ would then be acceptable, while one rising significantly faster would imply 
a decreasing interaction cross section o"int(p±o)) and by implication a decreasing cxtot; in 
contradiction with data. The DL fit to atot [7] gives e ^ 0.08, which would imply (at most) 
a p±o dependence like s° °^ = E^^^. However, cxtot already represents the unitarization of 
multiple-pomeron exchanges, and the 'bare' pomeron intercept should be larger than this, 
exactly by how much being a matter of some debate [65]. A value like ebarc ~ 0.12 is here 
near the lower end of the sensible range; the xg{x) shape is consistent with a rather larger 
value. Since it is actually the bare pomeron that corresponds to a single interaction, an 
-^cm^ behaviour is thereby acceptable. 
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3 Impact- Parameter Dependence 



In the simplest multiple-interactions scenarios, it is assumed that the initial state is the 
same for all hadron collisions. More realistically, one should include the possibility that 
each collision also could be characterized by a varying impact parameter b [1]. Within the 
classical framework we use here, b is to be thought of as a distance of closest approach, not 
as the Fourier transform of the momentum transfer. A small b value corresponds to a large 
overlap between the two colliding hadrons, and hence an enhanced probability for multiple 
interactions. A large b, on the other hand, corresponds to a grazing collision, with a large 
probability that no parton-parton interactions at all take place. 

In order to quantify the concept of hadronic matter overlap, one may assume a spher- 
ically symmetric distribution of matter inside a hadron at rest, p(x) d^x = p{r) d^x. For 
simplicity, the same spatial distribution is taken to apply for all parton species and mo- 
menta. Several different matter distributions have been tried. A Gaussian ansatz makes 
the subsequent calculations especially transparent, but there is no reason why this should 
be the correct form. Indeed, it appears to lead to a somewhat too narrow multiplicity 
distribution and too little of a pedestal effect. The next simplest choice, that does provide 
more fluctuations, is a double Gaussian 

p(r) oc — 3— exp <^ — 2 r + ~ 1 — 2 r • (12) 

This corresponds to a distribution with a small core region, of radius 02 and containing a 
fraction f3 of the total hadronic matter, embedded in a larger hadron of radius ai. If we 
want to give a deeper meaning to this ansatz, beyond it containing two more adjustable 
parameters, we could imagine it as an intermediate step towards a hadron with three disjoint 
core regions ('hot spots'), reflecting the presence of three valence quarks, together carrying 
the fraction /3 of the proton momentum. One could alternatively imagine a hard hadronic 
core surrounded by a pion cloud. Such details would affect e.g. the predictions for the t 
distribution in elastic scattering, but are not of any consequence for the current topics. 

For a collision with impact parameter b, the time-integrated overlap 0{b) between the 
matter distributions of the colliding hadrons is given by 



(1-/3)2 r 62 2p{i-p) r b^ \ 



The necessity to use boosted p(x) distributions has been circumvented by a suitable scale 
transformation of the z and t coordinates. The overlap function 0{b) is closely related to 
the Q{b) of eikonal models [24], but is somewhat simpler in spirit. 

The larger the overlap 0{b) is, the more likely it is to have interactions between partons 
in the two colliding hadrons. In fact, to first approximation, there should be a linear 
relationship 

(n(6)) = kO{b) , (14) 

where n = 0, 1, 2, . . . counts the number of interactions when two hadrons pass each other 
with an impact parameter b. At this stage k is an undefined constant of proportionality, to 
be specified below. 
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For each given impact parameter, the number of interactions n is assumed to be dis- 
tributed according to a Poissonian, Va = (n)^ exp{—{n))/n\, before energy-momentum 
and other constraints are included. If the matter distribution has a tail to infinity (as the 
Gaussians do), events may be obtained with arbitrarily large b values. In order to obtain 
finite total cross sections, it is necessary to assume that each event contains at least one 
semi-hard interaction. (Unlike the simpler, impact-parameter-independent approach above, 
where p± = no-interaction events are allowed as a separate class.) The probability that 
two hadrons, passing each other with an impact parameter b, will produce a real event is 
then given by 

oo 

VUb) = = 1 - = 1 - exp(-(n(6))) = 1 - exp{-kOib)) , (15) 

n=l 

according to Poisson statistics. The average number of interactions per event at impact 
parameter b is now {n{b)) = {n{b))/Vint{b), where the denominator comes from the removal 
of hadron pairs that pass without interaction, i.e. which do not produce any events. While 
the removal of n = from the potential event sample gives a narrower-than-Poisson in- 
teraction distribution at each fixed b, the variation of {n{b)) with b gives a 6-integrated 
broader-than-Poisson interaction multiplicity distribution. 

Averaged over all b the relationship (n) = Uint/cnd should still hold. Here, as before, (Tint 
is the integrated interaction cross section for a given regularization prescription at small 
p±, while the inelastic nondiffractive cross section and is taken from parametrizations [7,10]. 
This relation can be used to solve for the proportionality factor k in eq. (14). Note that, 
since now each event has to have at least one interaction, (n) > 1, one must ensure that 
Cint > o"nd- The p±o parameter has to be chosen accordingly small — since now the concept 
of no-interaction \o'w-p± events is gone, aesthetically it is more appealing to use the smooth 
P±o turnoff than the sharp p_Lmin cutoff, and thereby populate interactions continuously all 
the way down to = 0. The whole approach can be questioned at low energies, since then 
very small p±o values would be required, so that many of the interactions would end up in 
the truly nonperturbative p± region. 

Technically, the combined selection of b and a set of scattering p±i values now becomes 
more complicated [1,9]. It can be reduced to a combined choice of b and p±i, according to 
a generalization of eq. (8) 



dV Oib) 1 da 

exp 



dp±id26 (C) and dp J 



Oib) /-^/^^ 1 da , , 
—dp^ 



{O) Jp, (Tnddp'j 



(16) 



The removal of the n = non-events leads to a somewhat special definition of the average [1] 

The subsequent interactions can be generated sequentially in falling pj_ as before, with the 
only difference that da/dp\ now is multiplied by 0{b)/{0), where b is fixed at the value 
chosen above. 

Note that this lengthy procedure, via p(r) and 0{b), is not strictly necessary: the prob- 
ability Vn for having n interactions could be chosen according to any desired distribution. 
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Figure 3: Overlap profile 0{b) for a few different choices. Somewhat arbitrarily the 
different parametrizations have been normalized to the same area and average b, i.e. 
same J 0{b) d^b and J bO{b) d^b. (Recall that we have not specified b in terms of any 
absolute units, so both a vertical and a horizontal scale factor have to be fixed for 
each distribution separately.) Insert shows the region 6 < 2 on a linear scale. 



However, with only Vn known and an n selected from this distribution, there is no obvious 
way to order the interactions in p_L during the generation stage, with lower-p_L interactions 
modified by the flavours, energies and momenta of higher-p^ ones. (This problem is partly 
addressed in ref. [38], by a post-facto ordering of interactions and a subsequent rejection of 
some of the generated interactions, but flavour issues are not easily solved that way.) 

There is also another issue, the parton- level pedestal effect, related to the transition 
from hard events to soft ones. To first approximation, the likelihood that an event contains 
a very hard interaction is proportional to nVn, since n interactions in an event means n 
chances for a hard one. If the requested hardest p± is gradually reduced, the bias towards 
large n dies away and turns into its opposite: for events with the hardest p± ^ the 
likelihood of further interactions vanishes. The interpolation between these two extremes 
can be covered if an impact parameter is chosen, and thereby an 0{b), such that one can 
calculate the probability of not having an interaction harder than the requested hardest 
one, i.e. the exponential in eq. (16). 

If the Gaussian matter distribution is the simplest possible choice, the double Gaussian 
in some respects is the next-simplest one. It does introduce two new parameters, however, 
where we might have preferred to start with only one, to see how far that goes. As an 
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alternative, we will here explore an exponential of a power of b 

0{b) (X exp {-b'^} (18) 

where d ^ 2 gives deviations from a Gaussian behaviour. We will use the shorthand 
ExpOfPow((i = . . .) for such distributions. Note that we do not present an ansatz for p{r) 
from which the 0{b) is derived: in the general case the convolution of two p is nontrivial. 
A peaking of 0{b) at 5 = is related to one of p(r) at r = 0, however. 

A lower d corresponds to an overlap distribution more spiked at 6 = and with a 
higher tail at large b, Fig. 3, i.e. leads to larger fluctuations. Specifically, the height of the 
6 = peak is related to the possibility of having fluctuations out to high multiplicities. 
To give some feeling, an exponential, ExpOfPow((i = 1), is not too dissimilar to the old 
Pythia double Gaussian, with (3 = 0.5 and 02/01 = 0.2. Conveniently, the Tune A 
double Gaussian, still with (3 = 0.5 but now 02/01 = 0.4, is well approximated in shape 
by an ExpOfPow((i = 1.35). Another alternative, commonly used, is to assume the matter 
distribution to coincide with the charge distribution, as gauged by the electric form factor 
Gsipl) = (1 +Pi/Ai^)"^ with = 0.71 GeV^. This gives an 0{b) oc {pbfK^{^ib), which 
also is close in form to Tune A, although somewhat less peaked at small b. 

As indicated above, there are two key consequences of a an overlap profile choice. One is 
the interaction multiplicity distribution and the other the parton-level pedestal effect. These 
two are illustrated in Figs. 4 and 5, respectively, for pp at 1.8 TeV, with p_Lo = 2.0 GeV as 
in Tune A. The three frames of each figure illustrate how momentum conservation effects 
suppress the probability to have an event with large multiplicity. This effect is even stronger 
now that each 

interaction is allowed to undergo full shower evolution, so that it carries away more of 
the available energy. In the figures, the default lower shower cut-off of 1 GeV has been used; 
obviously a larger cut-off would give results intermediate to the two lower frames. Further, 
the possibility of two hard-scattering partons being part of the same shower is not included. 
Note that the suppression of the high-multiplicity tail implies that a distribution with large 
fluctuations in reality will have fewer interactions on the average than a less-fluctuating one, 
if they (as here) start with the same assumed average before the momentum conservation 
effects are considered. This means that the choice of p±_o is somewhat dependent on the 
one of overlap profile. 

Let us now study the hadron-level multiplicity distribution, and begin with UA5 data 
at 200 and 900 GeV [48]. Tune A then does impressively well. Fig. 6, in spite of primarily 
having been tuned to pedestal effects rather than multiplicity distributions. In this com- 
parison, we do not put too much emphasis on the low-multiplicity end, which is largely 
probing diffractive physics. Here the Pythia description is known to be too simple, with 
one or two strings stretched at low and no hard interactions at all. More relevant is the 
mismatch in peak position, which mainly is related to the multiplicity in events with only 
one interaction. Assuming that most hadronization parameters are fixed by e"'"e~ data, it is 
not simple to tune this position. The beam remnant structure does offer some leeway, but 
actually the defaults are already set towards the end of the sensible range that produces 
the lower peak position, and still it comes out on the high side. 

However, the main impression is of a very good description of the fluctuations to higher 
multiplicities, better than obtained with the old parameters explored in [1]. Of course, many 
aspects have changed significantly since then, such as the shape of parton densities at small 
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Figure 4: Distribution of the number of interactions for different overlap profiles 0{b), 
for pp at 1.8 TeV, top without momentum conservation constraints, middle with such 
constraints included but without (initial-state) showers, and bottom also with shower 
effects included. 
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Figure 5: Average number of interactions as a function of the p± of the hardest inter- 
action, for pp at 1.8 TeV, top without momentum conservation constraints, middle 
with such constraints included but without (initial-state) showers, and bottom also 
with shower effects included. 
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Charged multiplicity distribution at 200 GeV 
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Charged multiplicity distribution at 900 GeV 
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Figure 6: Charged multiplicity distribution at 200 and 900 GeV; different overlap 
profiles compared with UA5 data [48]. 
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Charged multiplicity distribution at 1800 GeV 
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Figure 7: Charged multiplicity distribution at 1800 GeV; different overlap profiles 
compared with E735 data [49]. 

X. One main difference is that Tune A 90% of the time picks subsequent interactions to be of 
the gg — *■ gg type with colour flow chosen to minimize the string length. Since each further 
interaction thereby contributes less additional multiplicity, the mean number of interactions 
can be increased, and this obviates the need for the more extreme double-Gaussian default 
parameters. 

If, nevertheless, one should attempt to modify the Tune A parameters, deviating from 
its ExpOfPow((i = 1.35) near equivalent, it would be towards a smaller d, i.e. a slight 
enhancement of the tail towards high multiplicities. An example is shown in Fig. 6, with 
ExpOfPow((i = 1.2) and p±Q = 1.9 GeV (at 1800 GeV, with the Tune A energy rescaling). 

However, the nice picture is shattered if one instead considers the E735 data at 1800 GeV 
[49], Fig. 7. Tune A gives a way too small tail out to large multiplicities, and also the 
ExpOfPow((i = 1.2) falls below the data. One would need something like an exponential 
with a rather low p±o = 1.6 GeV to come near the E735 data, and that then disagrees with 
the lower-energy UA5 data. Fig. 6. The agreement could be improved, but not to the level 
of Tune A, by playing with the energy dependence of p±o- However, the E735 collaboration 
itself notes that results from the two collaborations are incompatible over the whole UA5 
energy range and especially at 546 GeV, where both have data [49]. Furthermore, we do 
not have the expertise to fully simulate E735 selection criteria, nor to assess the impact of 
the large acceptance corrections. E735 only covered the pseudorapidity range |?7| < 3.25, 
so about half of the multiplicity is obtained by extrapolation from the measured region for 
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the 1800 GeV data. UA5 extended further and observed 70%-80%, depending on energy, 
of its muhiphcity. 

Obviously new experimental studies would be required to resolve the UA5-E735 ambi- 
guity. As it stands, presumably a tune adjusted to fit E735 would give disagreement with 
the CDF data that went into Tune A. In this particular case we suspect the differences to 
be of an experimental origin, but in other cases it could well be that the Pythia model 
is incapable of fitting different (correct) distributions simultaneously, the model not being 
perfect. Indeed, speaking in general terms, that is a main reason why we try to improve the 
model in this article. In this particular case and for the moment being, however, we choose 
to use the UA5-compatible Tune A as a convenient reference for a realistic multiplicity 
distribution at Tevatron energies. 

4 Correlations in Momentum and Flavour 

Consider a hadron undergoing multiple interactions in a collision. Such an object should be 
described by multi-parton densities, giving the joint probability of simultaneously finding n 
partons with flavours /i, . . . , /„, carrying momentum fractions Xi, . . . ,Xn inside the hadron, 
when probed by interactions at scales Qf, . . . ,Ql^. However, just like the standard one- 
particle-inclusive parton densities, such distributions would involve nonperturbative initial 
conditions that ultimately would have to be pinned down by experiment. We are nowhere 
near such a situation: the experimental information on double parton scattering, n = 2, 
boils down to one single number, the cxefr of eq. (10), and for n > 3 there is no information 
whatsoever. Wishing to make maximal use of the existing (n = 1) information, we thus 
propose the following strategy. 

As described above, the interactions may be generated in an ordered sequence of falling 
p±. For the hardest interaction, all smaller p_i_ scales may be effectively integrated out of 
the (unknown) fully correlated distributions, leaving an object described by the standard 
one-parton distributions, by definition. For the second and subsequent interactions, again 
all lower-p_L scales can be integrated out, but the correlations with the first cannot, and so 
on. 

The general situation is depicted in Fig. 8. This illustrates how, for the i'th interaction, 
only the correlations with the i — 1 previous interactions need be taken into account, with 
all lower p_L scales integrated out. Note, however, that this is only strictly true for the 
hard scatterings themselves. The initial-state shower evolution of, say, the first interaction, 
should exhibit correlations with the z'th at scales smaller than p±i. Thus, the p± ordering 
(or equivalently, a virtuality ordering) is in some sense equivalent to a time ordering, with 
the harder physics being able to influence the softer physics, but not vice versa. For 
two interactions of comparable p± this order may appear quite arbitrary, and also should 
not matter much, but consider the case of one very hard and one very soft interaction. 
The soft one will then correspond to a long formation time (field regeneration time) [66], 
~ p/pi ~ l/p±, and indeed it is to be expected that the hard one can pre-empt or at least 
modify the soft one, whereas the influence in the other direction would be minor. This 
gives additional motivation to the choice of a p± ordering of interactions. 

The possibility of intertwined shower evolution is not (yet) addressed. Rather, we intro- 
duce modified parton densities, that correlate the z'th interaction and its shower evolution 
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Figure 8: Schematic representation of the evolution of parton shower initiators in a hadron 
coUision with n interactions (see text). 



to what happened in the 2 — 1 previous ones, but we do not let the previous showers be 
affected by what happens in subsequent interactions. As partons are successively removed 
from the hadron by hard scatterings at smaller and smaller p_i_ scales, the flavour, mo- 
mentum and colour structure of the remaining object changes. The colour structure in 
particular is a thorny issue and will be discussed separately, in the next Section. Here, we 
focus on deriving a set of parton distributions for a hadron after an arbitrary number of 
interactions have occurred, on the need for assigning a primordial transverse momentum to 
shower initiators, and on the kinematics of the partons residing in the final beam remnants. 

Our general strategy is thus to pick a succession of hard interactions and to associate 
each interaction with initial- and final-state shower activity, using the parton densities 
introduced below. The initial-state shower is constructed by backwards evolution [67], 
back to the shower initiators at some low Qq scale, the parton shower cutoff scale. Thus, 
even if the hard scattering does not involve a valence quark, say, the possibility exists that 
the shower will reconstruct back to one. This necessitates dealing with quite complicated 
beam remnant structures. For instance, if two valence quarks have been knocked out of 
the same baryon in different directions, there will be three quarks, widely separated in 
momentum space, of which no two may naturally be collapsed to form a diquark system. 

In the old model, technical limitations in the way the fragmentation was handled made it 
impossible to address such remnant systems. Consequently, it was not possible to associate 
initial-state radiation with the interactions after the first, i.e. the one with the highest p± 
scale, and only a very limited set of qq and gg scatterings were allowed. 

In a recent article [3], the Lund string model was augmented to include string systems 
carrying non-zero baryon number, by the introduction of 'junction fragmentation'. In the 
context of multiple interactions, this improvement means that almost arbitrarily compli- 
cated beam remnants may now be dealt with. Thus, a number of the restrictions that were 
present in the old model may now be lifted. 
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4.1 Parton Densities 

As mentioned above, we take the standard parton density functions as our starting point 
in constructing parton distributions for the remnant hadronic object after one or several 
interactions have occurred. Based on considerations of momentum and flavour conservation 
we then introduce successive modifications to these distributions. 

The first and most trivial observation is that each interaction i removes a momentum 
fraction Xi from the hadron remnant. This is the fraction carried by the initiator of the 
initial-state shower, at the shower cutoff scale Qq, so that the two initiators of an interaction 
together carry all the energy and momentum eventually carried by the hard scattering and 
initial-state shower products combined. To take into account that the total momentum of 
the remaining object is thereby reduced, already in the old model the parton densities were 
assumed to scale such that the point x = 1 would correspond to the remaining momentum 
in the beam hadron, rather than the total original beam momentum, cf. eq. (9). In addition 
to this simple x scaling ansatz we now introduce the possibility of genuine and non-trivial 
changes in both shape and normalization of the distributions. 

4.1.1 Valence Quarks 

Whenever a valence quark is knocked out of an incoming hadron, the number of remaining 
valence quarks of that species should be reduced accordingly. Thus, for a proton, the 
valence d distribution is completely removed if the valence d quark has been kicked out, 
whereas the valence u distribution is halved when one of the two is kicked out. In cases 
where the valence and sea u and d quark distributions are not separately provided from 
the PDF libraries, we assume that the sea is fiavour-antifiavour symmetric, so that one can 
write e.g. 



Here and in the following, [qs) denotes the q valence (sea) distribution. The parametrized 
u and u distributions should then be used to find the relative probability for a kicked-out 
u quark to be either valence or sea. Explicitly, the quark valence distribution of flavour 
/ after n interactions, qf^n{x,Q'^), is given in terms of the initial distribution, Qf^oi^^Q'^), 
and the ratio of remaining to original qj valence quarks, Nf^n/ Nf^Q, as: 



where N^vo = 2 and A^^vo = 1 for the proton, and x G [0, X] is the fraction of the original 
beam momentum {^2^=1 is the total momentum fraction already taken out of the incoming 
hadrons by the preceding parton-shower initiators). The dependence of gjvn is inherited 
from the standard parton densities Qf^o, and this dependence is reflected both in the choice 
of a hard scattering and in the backwards evolution. The factor 1/X arises since we squeeze 
the distribution in x while maintaining its area equal to the number of q/ valence quarks 
originally in the hadron, A^/vo, thereby ensuring that the sum rule. 
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is respected. There is also the total momentum sum rule, 




qfn{x, Q"^) + gn{x, Q"^) \xdx = X. 



(22) 



Without any further change, this sum rule would not be respected since, by removing a 
valence quark from the parton distributions in the above manner, we also remove a total 
amount of momentum corresponding to (x/v), the average momentum fraction carried by 
a valence quark of flavour /: 



has already been taken into account by the 'squeezing' in x of the parton distributions (and 
expressed in eq. (22) by the RHS being equal to X rather than 1). By scaling down the 
qv distribution, we are removing an additional fraction, (xf^n), which must be put back 
somewhere, in order to maintain the validity of eq. (22). 

Strictly speaking, (x/vo) of course depends on which specific PDF set is used. Never- 
theless, for the purpose at hand this variation is negligible between most modern PDF sets. 
Hence we make the arbitrary choice of restricting our attention to the values obtained with 
the CTEQ5L PDF set [6]. 

More importantly, all the above parton densities depend on the factorization scale Q^. 
This dependence of course carries over to (x/vo), for which we assume the functional form 



inspired by the ds = dlog(log(5^/A^) oc dQ'^/Q'^ dsiQ^) pace of evolution, where d{x)/ds ~ 
—B{x) suggests a solution of the form (x) oc exp{—Bs) ^1/ (l + Bs). Reasonable fits to the 
CTEQ5L valence quark distributions in the proton are obtained for = 0.385, = 1.60, 
= 0.48 and B^ = 1.56, with the isospin conjugate for neutrons. 

Essentially nothing is known about parton densities for other baryons, such as the 
reasonably long-lived hyperons A°, S"*"'^, and , which can undergo secondary inter- 
actions that one may wish to study. We here use essentially the same parton densities and 
parameters as for protons. Thus the influence of the larger strange quark mass is neglected, 
which ought to lead to harder x spectra for s quarks and softer for everything else. The 
fact that the two proton valence u quarks have a harder distribution than the single d one 
is carried over to other baryons with two equal quarks, while the average (sum) of the u 
and d distributions are used for baryons with three unequal (equal) quarks. For mesons 
and the Vector Meson Dominance part of photons one could use a similar strategy, with 
vr"*" measurements as a starting point instead of protons, while the anomalous part of the 
photon densities is perturbatively calculable. There are still many further assumptions that 
would have to go into a complete model of multiple interactions in 7p and 77 events [12], 
however, and so far we did not pursue this further. 

We now know how much momentum is 'missing' in eq. (22). It is not possible to put 
this momentum back onto valence quarks without changing the shape of the distributions 
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1 + Bf log (log(max(Q2, 1 GeV')/Alc^)) ' 
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(beyond the mere 'a; squeezing') or invalidating eq. (21). Rather, we here assume that the 
missing momentum is taken up by the sea+gluon distributions, which thus are scaled up 
slightly when a valence quark is kicked out. This enhancement of the sea+gluon momentum 
fraction may over- or undercompensate the 'x squeezing' reduction, depending on whether 
the kicked-out valence quark had a small or large x. However, before the procedure can be 
discussed in more detail, we must consider another effect which affects the normalization 
of the sea: changes in the content of the sea itself. 



4.1.2 Sea Quarks and their Companions 

When a sea quark is kicked out of a hadron, it must leave behind a corresponding antisea 
parton in the beam remnant, by flavour conservation. We call this a companion quark. In 
the perturbative approximation the sea quark qg and its companion q^ come from a gluon 
branching g — > qs + qc, where it is implicitly understood that if qs is a quark, qc is its 
antiquark, and vice versa. This branching often would not be in the perturbative regime, 
but we choose to make a perturbative ansatz, and also to neglect subsequent perturbative 
evolution of the qc distribution. Even if approximate, this procedure should catch the key 
feature that a sea quark and its companion should not be expected too far apart in x (or, 
better, in Inx). 

With this approximation, we obtain the distribution from the probability that a sea 
quark qg, carrying a momentum fraction Xg, is produced by the branching of a gluon with 
momentum fraction y, so that the companion has a momentum fraction x = y — Xs, 

qdx; X,) = C g{y) Pg^q,qe(2) 5(xs - zy) dz 







C 9{y) Pg. 



X«\ 1 



y J y 



^ C?i^p(^\ (25) 

Xs + X \Xs + xJ 

with C a normalization constant to be determined below, and Pg^q^q^ the DGLAP splitting 
kernel 

Pg^qeqcW = I{z'+{1-Zf). (26) 

In view of the approximate nature of the procedure, allowing a generic g{x) shape would 
give disproportionately complex expressions. Instead, the following simple ansatz for the 
gluon distribution at low is used: 

(1 — x)P 

9{x) oc ^ ^, (27) 

X 

with the integer choices p = 0, 1, 2, 3, 4 giving a range of variability for the large-x behaviour 
of the distribution and the 1/x controlling the small-x behaviour. Note that all the above 
equations are defined assuming no previous energy loss and, as for the valence quarks, 
should be 'squeezed' by a factor X, to ensure momentum conservation. 
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Companion Distributions 




Figure 9: Companion distributions for p = 4 (solid lines) and p = (dotted lines), for two 



different values of Xg. 



The overall normalization of a companion quark distribution is obtained by imposing 
the sum rule: 



l-Xs 



qc{x;Xs) dx = 1. (28) 
Inserting eqs. (25)-(27) and inverting, one obtains the normalization constants Cp-. 



Co = 


3xs 


(29) 


2 3/g(3 ~1~ 2xg) 


Ci = 


3xs 


(30) 


2 - x2(3 - Xg) + 3xg log(xs) 


C2 = 


Sxg 


(31) 


2(1 - Xg)(l + 4x, + 4x2) + 6xg(l + Xg) log(xg) ' 


C3 = 


6xg 


(32) 


4 + Xg(27 - 31x2) + 6xg(3 + 2xg(3 + Xg)) log(xg) ' 


c, = 


3xg 


(33) 


2(1 + 2xg)((l - Xg)(l + Xg(10 + Xg)) + 6xg(l + Xg) log(xg)) ' 



To illustrate. Fig. 9 shows properly normalized companion momentum distributions for 
p = 4 and p = 0, each for two different values of Xg. There are two noteworthy aspects 
about these distributions. Firstly, for p = the discontinuity at the point x = 1 — Xg is 
merely an artifact of our parametrization of the gluon density, i.e. that g{x) does not vanish 
at X = 1. That problem is absent for more realistic p values; in our continued discussions we 
will use p = 4 as default, this being the closest to the CTEQ5L small-Q^ gluon distribution. 
Secondly, the falling gluon distribution convoluted with the almost flat g ^ qq splitting 
kernel give distributions that roughly tend to a constant qc{x;Xs) ~ Cp/2x^ below Xg and 
exhibit power-like fall-offs qc{x;Xs) oc l/x2 above it, with some modulation of the latter 
depending on p. In order to display the probability per logx interval. Fig. 9 gives xQc 
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rather than q^. itself, and then a peaking occurs around should be expected from 

the symmetric splitting kernel in eq. (26). 

Also here, the question arises of ensuring that the total momentum sum rule, eq. (22), 
is respected, but now the difference has the opposite sign; by adding a companion quark 
distribution, we are in effect bookkeeping a part of the flavour and momentum content of 
the sea separately. One possibility is that this momentum comes only from the sea+gluons 
and that the valence quarks are not affected, i.e. that the rest of the sea+gluons fluctuate 
down, in order to compensate. 

The amount of momentum that will have to be compensated for each companion quark, 
{xcn) = ^{xco), with (xcn) defined in analogy with eq. (23), is straightforward to compute 
using the distribution in eq. (25) and the normalizations given by eqs. (29)-(33): 

_ 5-9a;s + 6x^ -2x^ + 3 log , , 

^'^^°^^=° - (x3-l)(2-x, + 2x2) ' ^"^^^ 

{xco}p=i = -1 - 3xs + - — — — , (35) 

2 — 3xf + + Sx^ log x^ 

. . _ 19 + 24xs - 39x,^ - + 6(1 + Gx^ + 4x^) logXs 

{Xco}p=2 - J -1 - 3x, + 3x2 + ^3 _ 3^^(i + i^g^.^ ' ^"^^^ 

-7 - 21xs + 15x^ + 13x^ - 2(1 + 9xs + Ux^ + 2x^) logXs 
(Xco)p=3 - 3xs 4 + 27x,-31x3 + 6x,(3 + 6xs + 2x2)logXs ' ^ ' 

3(5 + 24xs - - 24:xl - x^) + 4(1 + 12xs + 2Axl + 8x^) logx^ 
(Xco)p=4 - 3xs ^ 2x,)i-l - 9xs + 9x2 + ^3 _ 6^^(i + log^^) • ^^^) 

4.1.3 Sea Quark and Gluon Density Normalizations 

As described above, the normalization of valence and companion distributions is fixed by 
the respective number of quarks, i.e. the sum rules (for each flavour, /) 

/ qfvnix,Q^) dx = Nf^n, (39) 
Jo 



/ qfcjn{x;xsj) dx = 1 (for each j) 

^0 



(40) 



where X is still the longitudinal momentum fraction left after the n previous interactions 
and A^/vn is the number of q/ valence quarks remaining. The index j on the companion 
distribution, qjc^n, counts different companion quarks of the same flavour, /. 

On the other hand, the sea+gluon distributions do not have fixed multiplicities, hence 
no corresponding sum rules exist for their normalizations. We use this freedom to fulfill the 
last remaining sum rule, eq. (22), letting the sea+gluon normalizations fluctuate up when 
we reduce a valence distribution and down when we add a companion distribution. In 
addition, the requirement of a physical x range is of course maintained by still 'squeezing' 
all distributions into the interval x G [0,X]. 

For simplicity, and since eq. (22) only furnishes us with one equation of constraint, we 
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assume the same scale factor for all sea flavours as well as for the gluon, i.e. 

qfsix,Q'^) ^ ag/s(x,(5^), (41) 
g{x,Q') ^ ag{x,Q^). (42) 

The momentum sum rule now reads: 



/ L j 

^^^" Qfvo{x, Q^) + X qfc,o{x; Xj) + ag/so(a;, Q^) 



N 



fvO 



agn{x, Q^) j X dx 
+ agQ{x, Q^) I X dx 



fvn 



N-. 



fvO 



Qfvoix, Q^) + qjcjoix; xj) 



X dx 



/ / / f,j 



and hence, 



(43) 



(44) 



1 - E/^/vo(a;/vo) 

One easily checks that a = 1 before the first interaction, as it should be, and that a is 
driven larger by 7V/vn < A''/vo, while introducing companion quarks drives it the opposite 
way, also as expected. 



4.2 Beam Remnants 

The longitudinal momenta and flavours of the initiator partons are defined by the sequence 
of p_L-ordered hard scatterings and their associated initial-state showers, as described above. 
What is left in the beam remnant is then a number of partons, with flavours given by 
the remaining valence content plus the number of sea quarks required for overall flavour 
conservation. That is, gluons in the remnant are not explicitly accounted for, but are 
imphcit as confinement clouds around the quarks and as unresolved originators of sea quark 
pairs. 

A remnant may thus contain several objects but, when the colour configuration is stud- 
ied, simplifications can occur. A colour antitriplet qq pair in the remnant can be associated 
with a diquark, a colour singlet qqq triplet with a baryon, and a colour singlet qq pair with 
a meson, see Fig. 10. When hadrons are formed, the standard string fragmentation relative 
probabilities are used to select spin and other quantum numbers, i.e. whether vr or p, etc. 

It is here assumed that the respective pair/triplet has a sufficiently small invariant mass 
that it can reasonably be projected onto a single composite state. Thus a qq system with 
large invariant mass would define a string that could fragment into several mesons, rather 
than collapse to a single meson. In principle this could be modeled dynamically, but it would 
require the introduction of some nonperturbative parameters, to describe the partitioning 
of the proton into arbitrary-mass subsystems. At this point, we consider it meaningful only 
to study a few specific scenarios for which partons to allow in the formation of composite 
objects. We have chosen four such: 
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c) 

Figure 10: Examples of the formation of composite objects in a baryon beam remnant: (a) 
diquark, [h) baryon and (c) meson. 

1. No composite objects are formed ah initio. All partons act as single units, either as 
endpoints (quarks) or kinks (gluons) on strings that fragment in the normal way. 

2. Composite objects may be formed, but only when all partons involved in the formation 
are valence quarks. 

3. The formation of diquarks may involve both valence and sea quarks, but the formation 
of colour singlet subsystems (i.e. hadrons) is still restricted to involve valence quarks 
only. 

4. Sea quarks may also be used for colour singlet formation. 

The idea is thus that (spectator) valence quarks tend to have comparable velocities, while 
sea quarks can be more spread out and therefore are less likely to form low-mass systems. 

Whether composite systems in the beam remnant are formed or not has important 
consequences for the baryon number flow. For pp collisions at 1.8 TeV CM energy, we show 
in Fig. 11 the Feynman x (left plot) and rapidity (right plot) distributions for the baryon 
which 'inherits' the beam baryon number. We denote this baryon the 'junction baryon'. To 
better illustrate what happens to each of the two initial beam baryon numbers separately, 
only distributions for the junction baryon, not antzbaryon, are shown. Possibilities 1 and 2 
above are compared with the old multiple interactions model (Tune A). One immediately 
observes that the beam baryon number migrates in a radically different way when diquark 
formation is allowed or not (compare the dashed and dotted sets of curves). In fact, in the 
new model it is not possible to reproduce the old distribution (compare the solid curve). 
This comes about since, even when all possible diquark formation is allowed in the new 
model, it is not certain that the beam remnant actually contains the necessary quark 
content, hence in some fraction of the events the formation of a beam remnant diquark is 
simply not possible. Here is thus an example where the introduction of more physics into 
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Figure 11: Feynman x (left) and rapidity (right) distributions for junction baryons: distri- 
butions are shown for Tune A of the old MI scenario (solid lines), and for the new model 
with diquark formation in the beam remnant switched on (dashed) and off (dotted). 



the model has given rise to a qualitatively different expectation: the beam baryon number 
appears to be stopped to a larger extent than would previously have been expected. 

One should note that, also at later stages, a small-mass string piece can collapse to a 
single hadron, as part of the normal string fragmentation procedure. There, however, it 
is intended only to cover a rare low-mass tail of systems mainly defined by hard processes 
and perturbative shower evolution, while the simplifications considered for the above kine- 
matical configurations are quite common and in a nonperturbative context. Beam remnant 
quarks that were not collapsed in the nonperturbative first stage could in the later stage 
be collapsed with other partons, e.g. from the showers. Also collapses not allowed in the 
more restrictive scenarios above could occur at this stage. Whether that happens or not 
depends on the transverse and longitudinal momenta that will be defined below. 

4.3 Primordial k±_ 

Until now, we have considered only the longitudinal part of parton momenta. In reality, 
partons are also expected to have some non-zero /c^ values caused by Fermi motion inside 
the incoming hadrons. This kind of k_[_ is denoted 'primordial /c^', since it is not generated 
by the (DGLAP) shower evolution nor from hard interactions, but rather represents an 
input to the perturbative stages of the event. Based on Fermi motion alone, one would 
expect values of the order of a few hundred MeV, just like in eq. (6). But to reproduce 
e.g. the p_[_ distributions of Z bosons produced in hadron-hadron collisions, one notes a need 
for a significantly larger nonperturbative input, either in parton showers or in resummation 
descriptions [68]. This problem is still awaiting a satisfactory explanation, although some 
ideas have been explored that might alleviate it [69]. Until such an explanation has been 
found, we therefore have reason to consider an effective 'primordial k±\ at the level of the 
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initiators, larger than the one above. For simphcity, a parametrized Q-dependent width 



a{Q) = max ^min, '^°^Y+~QTJQ j 

is introduced, where a is the width of the two-dimensional Gaussian distribution of the 
initiator primordial k± (so that {k"^) = cr^), Q is the scale of the hard interaction, a^o is 
the value asymptotically approached as Q ^ oo, and Qi is the Q scale where a = t^cToc- A 
reasonable fit to the few available experimental 'data points', 

o-{Q ~ 1 GeV) ~ cTjnin ~ 0.36 GeV (fragmentation) 
a(g~5-10GeV) ^ 0.9 GeV (EMC) [70] (46) 
a{Qr^Mz) ^ 2 GeV (Tevatron) [68], 

is obtained with the values aoo = 2.1 GeV and Qi = 7 GeV. The amin in eq. (45) represents 
a minimum broadening, at the level of Fermi motion, which we take to be the standard 
fragmentation p± width. In addition to partons participating in relatively soft interactions, 
this minimum broadening is also applicable to the remnant partons, which by definition do 
not participate in hard interactions and hence are not naturally associated with a particular 
Q scale. 

Apart from the selection of each individual kj_, there is also the requirement that the 
total k± of the beam adds up to zero. The question of how partons recoil off one another in 
transverse momentum space inside hadrons is so far largely unaddressed in the literature. 
We imagine a few different possibilities here: 

1. The primordial transverse momenta are generated at a stage where the partons have 
low virtualities and hence large wavefunctions. Moreover, at least for that part which 
is due to Fermi motion, the dynamics responsible for the generation of primordial k± 
is that of a Fermi gas of partons in equilibrium, where each parton has received its 
total primordial k_i_ through a sequence of many collisions with many different partons. 
Therefore, one possibility is to let the recoil of one parton be shared uniformly among 
all other initiator and remnant partons. 

2. Since Fermi motion alone appears unable to account for the bulk of primordial k± 
in large-Q interactions, there may exist a mechanism, involving presently unknown 
dynamics, which ties the generation of this k± to the presence of a large virtuality 
in the interaction, for instance by unresolved/unresummed bremsstrahlung radiation 
off the initiator parton and/or by the k± compensations in the shower being of a 
not strictly local nature (as happens e.g. in the dipole description of parton showers 
[36]). Regardless of the exact nature of the mechanism, recoils should in this picture 
primarily be taken up by initiators and beam remnant partons which are close in 
colour space. 

(a) The extreme variant is here to let the recoil of a particular parton be taken up by 
its nearest colour neighbours only. 

(b) A perhaps more realistic possibility is to let the compensation happen along a 
parton chain in colour space, with successive dampening of the compensation along 
the chain. 
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These three possibihties are included in the present study. However, the possibihties involv- 
ing colour chains are complicated by the fact that the colour connections between initiator 
and remnant partons are very poorly known, since no perturbative information is available. 
These problems will be discussed in more detail in Section 5 below. 

Irrespective of which particular method is used to ensure k±i = 0, the question 
now arises how the kinematics of the initially collinear partons should be reinterpreted to 
include non-zero k± assignments. This may be done either by associating the generation 
of k± with the building up of space-like virtualities among the partons, or by keeping 
the partons massless while allowing (non-perturbatively small) longitudinal momentum 
transfers between the beam remnants. (In the latter approximation also heavy quarks 
are kept massless rather than assigned a spacelike virtuality; when initial-state showers 
are included no heavy quarks need be assigned to the beam remnants, however.) In the 
first case, the invariant mass of initiator and beam remnant partons combined in each 
hadron is maintained equal to the original hadron mass, while in the second the mass can 
be significantly larger. The difference, however, should be considered mostly technical, 
since the momentum transfers involved are quantitatively small. For this study, the second 
option is chosen, since this avoids potential technical problems in dealing with string systems 
having negative mass squares. 

For a specific interaction, consider a pair of massless initiator partons in their rest frame, 
before k± is added: 



if the system should still be at longitudinal rest. Since we are merely reinterpreting the 




(47) 



With primordial k± included, these momentum vectors should now be recast as 




(48) 



kinematics of the initial-state partons, the centre-of-mass en 
should be left unchanged. To ensure this, it is simple to solve 





(49) 



to obtain the required clS cl function of p±2, and s: 




(50) 



with A the standard Kallen function 




(51) 



Naturally, only k± assignments which result in > are acceptable. 
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4.4 Beam Remnant Longitudinal Momenta 

In addition to flavours and transverse momenta, the beam remnants must also together 
carry the remaining fraction, approximately X (as defined by eq. (20)) of longitudinal 
momentum. The sharing is based on the character of the remnant constituents. First a 
fraction x is defined for each constituent, and then these x fractions are rescaled for overall 
energy and momentum conservation. 

Thus a valence quark receives an x picked at random according to a small-Q^ valence- 
like parton density, proportional to (1 — x)"" j \fx^ where a = 2 for a u quark in a proton 
and a = 3.5 for a d quark. A sea quark must be the companion of one of the initiator 
quarks, and can have an x picked according to the qdx^x^ distribution introduced above. 
In the rare case that no valence quarks remain and no sea quarks need be added for flavour 
conservation, the beam remnant is represented by a gluon, carrying all of the beam remnant 
longitudinal momentum. 

Among composite objects, a diquark would naively obtain an x given by the sum of its 
constituent quarks, while baryons and mesons would receive an x equated with the z value 
obtainable from a fragmentation function, in this study the Lund symmetric fragmentation 
function. However, earlier studies on quark-diquark remnants [1] have shown that, within 
the multiple interactions formalism, it is very difficult to accommodate observed remnant 
multiplicity distributions if the composite system (the diquark) does not take a much larger 
fraction than implied by the naive estimate above. Physically, this could correspond to the 
momentum carried by a surrounding pion/gluon cloud being larger for a composite object 
than for a single parton. The possibility of enhancing the x values picked for composite 
objects is therefore retained in the present study. 

Finally, once x values (and primordial /c_l) have been picked for each of the remnants, an 
overall rescaling is performed such that the remnants together carry the desired longitudinal 
momentum. The simplest way to accomplish this would be to fix the normalization of the 
beam remnant x values on each side separately, by requiring conservation of longitudinal 
momentum, Xj^j^ = X on each side. Unfortunately, the introduction of non-zero 
values with massless partons and on-shell hadrons rules out such a simple approach, since 
energy would then not be conserved. Instead, small non-zero lightcone momentum fractions 
in the direction opposite to the parent hadron direction must be allowed. As already noted 
in Section 4.3, this procedure should be thought of merely as a technical trick, necessitated 
by insisting on a description in terms of on-shell partons. 

The amount of light-cone momentum removed from the remnant system by each pair of 
initiators, i, is in the overall cm frame of the event (omitting subscript i to avoid cluttering 
the notation), 

= E-^^vT = lil + P.){E[ + E',) , 

where subscripts 1 and 2 denote the initiator partons of the hard scattering on each side 
respectively, and the primed frame is chosen as the longitudinal rest frame, defined by 
eq. (48). Then the boost is only along the z direction, 

P.=P = ^^, (53) 

Xi + X2 
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and from equations (48)-(50) 

E[ + E',= 

Inserting these results in eq. (52) above yields 



(54) 



w~ 




(55) 



For vanishing p± this simplifies to the familiar s± = s = X1X2S, = Xiy/s and w = X2\fs- 
The light-cone momenta remaining for the combined beam remnant system are thus: 



w~ 



i 



s — 




W~ . 

rem rem 



(56) 

(57) 
(58) 



In extreme cases, it may happen that the hard interactions have removed so much en- 
ergy and momentum from the beam remnants that the remnant system nominally becomes 
space-like, VF^em < O5 if large k± values have been assigned. Though not strictly speaking 
unphysical, such a situation could lead to problems at the fragmentation stage. The re- 
quirement W^^^ > is therefore imposed as an additional constraint when primordial k± 
values are assigned. 

The X values picked for the beam remnant partons are now interpreted as fractions of 
the light-cone momenta, W^^^^ and W~j^, of the beam remnant system, modulo an overall 
rescaling on each side, to leave room for overall momentum conservation. Using index j to 
refer to beam remnant partons on side 1 and index k for the ones on side 2, we thus make 
the identification 
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(59) 
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where 
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-p\ and a and (3 are global normalizations, to be determined from overall 



energy and momentum conservation in the beam remnant system: 

E^/ + E^^ = "^rtmE^^- + T^E 
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Equating these expressions with eqns. (56) and (57), one obtains for a and /3 

" = ^vi::T-^, ■ 

n rem ^'^'^2 yv i ^ y\ \yy rem? "^"^ 1 ; "^"^ 2 / /'«Q^ 

with the Kallen function given by eq. (51) and Wi and W2 the total transverse masses of 
each of the two beam remnant systems, 

Finally, for physical choices of Xj^k and primordial /c^, the sum of the individual remnant 
transverse masses must be smaller than that of the total remnant system, Wi + W2 < H^rem- 
If this is not the case, new k± sets and/or new xj^k values are tried, until a physical set of 
values is found. 




5 Colour Correlations and String Topologies 

The formalism described in the previous Sections may be used to obtain a sequence of 
hard scatterings with associated initial- and final-state showers and ordered by the p± 
of the hard interactions they contain. Kinematics is completely specified for all partons 
involved in the scatterings, in the associated showers as well as in the left-behind beam 
remnants. However, at some point, the time evolution of this system results in inter-parton 
distance scales larger than about a femtometer, where the perturbative QCD description 
in terms of partons breaks down and must be supplanted by one of hadrons. Now, if each 
parton hadronized independently of the rest of the event, the information on kinematics 
and flavours alone would suffice to pass from the language of partons to that of hadrons, 
but confinement is precisely not a strictly local effect. Rather, it is a statement about a 
(colour singlet) system of partons, and hadronization is one concerning the evolution when 
colour charges inside such a system have been imparted with large momenta relative to 
each other. It is therefore not meaningful to study the hadronization of a single parton in 
isolation. Instead, it becomes necessary to consider the interplay between colour charges 
and to take correlations into account when modeling the hadronization process. 

In this Section, we begin by considering the hadronization of systems containing non- 
zero baryon numbers, but where the colours of all participating partons are known. We 
will use the planar approximation of QCD as a starting point, with a junction picture for 
the baryons. 

Next, since the hard-scattering and parton-shower histories discussed above do not 
provide sufficient colour information — specifically information about how the different 
scattering initiators are correlated in colour is completely absent — we consider several 
possibilities for the assignment of correlated colours to the parton-shower initiators of the 
scatterings. Combined with the hadronization model, this allows us to study what is 
obtained under 'minimal' assumptions. 
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Finally, further issues are whether the original colour arrangement survives all the way to 
the long-distance hadronization era, and whether nonlinear effects arise in the hadronization 
process itself. That is, with several partons and string pieces moving out from the collision 
process, these partons and pieces will largely overlap in space and time. We do not know 
whether such overlaps can lead to colour rearrangements or nontrivial hadronization effects, 
e.g. of the Bose-Einstein kind. In principle e"'"e" — > W"*"W^ — ^ ^1^2*13^4 offers a clean 
environment to study such crosstalk, but experimental results are inconclusive [79]. In 
hadronic collisions, Bose-Einstein studies by UAl and E735 also give a splintered image [80]: 
the strength parameter of BE effects drops with increasing particle density, consistent with 
a picture where a higher multiplicity comes from having several independently hadronizing 
strings, but the BE radius also increases, which suggests correlations between the strings. 



5.1 Hadronization 

Taking colour interference into account when modeling the hadronization process could 
easily become an unmanageable task. One simplification (disregarding baryons for the 
moment) is to go to the limit with infinitely many colours, Nc oo [71]. In this limit 
the confinement force acting on a gluon is twice that on a quark, i.e. the gluon colour 
and anticolour charges decouple. Further, colour diagrams are planar, so that final-state 
colour-anticolour pairs are always uniquely matched, via an unbroken colour line through 
the diagram. 

In the Lund string model [4], the two ends of such a colour line define a string piece. 
The string piece can be viewed as a Lorentz covariant and causal implementation of a linear 
confinement potential between the two partons. Transverse degrees of freedom here play no 
dynamical role, but one can visualize the colour field lines as compressed into a tube of a 
typical hadronic width (~ 1 fm). As the partons move apart and a string piece is stretched 
out, it can break by the production of new qq pairs that screen the endpoint colour charges. 
The q of one such break and the q of an adjacent break together define a meson, which 
may be unstable and decay further. 

The classical example is e^e~ — > qqg, where one may assign a red colour to the q, 
antired-fgreen to the g and antigreen to the q, so that the string consists of two pieces, one 
q-g and one g-q. There is no piece directly between q and q, with observable consequences 
in the particle flow [72]. 

Turning now to baryon beams as the more interesting and difficult example, we picture 
the initial state of a baryon as consisting of three valence quarks connected in colour via a 
central 'junction', cf. Fig. 12. At the most basic level, such a picture finds its motivation by 
considering the simplest locally gauge invariant operator in SU (3) which carries non-zero 
baryon number [73]: 



,-1 L V JV(x,xi) J 



(65) 



Physically, this operator assigns the space-time coordinates Xi to three valence quarks (with 
fiavours and colours a,) and connects each of them via the gluon field along the path V 
to the point x (with P representing the path-ordering operation), which we may identify as 
the locus of the string junction. Such ideas were already introduced in the early string model 
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• Valence quark in beam hadron 



Figure 12: The initial state of a baryon, consisting of 3 valence quarks connected antisym- 
metrically in colour via a central 'string junction', J. 



of hadrons [73-75], and have been used to construct baryon wavefunctions in confinement 
studies [76]. In a recent article [3], we argued that this picture also arises naturally from 
string energy minimization considerations. 

In a collision, the fact that a gluon carries colour implies that the junction will, in 
general, be separated in colour space from the original valence quarks. As a simple example, 
consider a valence qq — > qq scattering in a pp collision. The exchanged gluon will flip 
colours, so that each junction becomes attached to a q from the other proton. 

Hence the junction may well end up colour-connected to partons — or chains of partons 
— which are widely separated in momentum space and of which no two may be naturally 
considered to form a diquark system. To describe the hadronization of such systems, a 
model capable of addressing colour topologies containing explicit non-zero baryon numbers, 
here in the form of junctions, becomes necessary. Such a model was first developed in [3], for 
dealing with the colour topologies that arise in baryon number violating supersymmetric 
scenarios. In the following, we show how this approach may be applied in a multiple 
interactions context to describe the physics of beam remnants. 

We begin by considering a simplified situation where only one of the initial beam par- 
ticles is a baryon. Leaving the ambiguities in assigning correlated colours aside for the 
moment. Fig. 13 gives an example of how the colour structure of a 7p collision might look. 
In this example, the final-state colour-singlet system containing the junction consists of the 
three string pieces J — q^, J — qvs, and J — g — q'. The two other string systems in the event, 
qvi — q' and qv2 — ^q^ are standard and do not concern us here. 

To understand how the junction system hadronizes, the motion of the junction must 
first be established. This can be inferred from noting that the opening angle between any 
pair of the connected string pieces is 120° in the rest frame of the junction, i.e. in that frame 
the system consisting of the junction and its nearest colour-connected neighbours looks like 
a perfect Mercedes topology. This is derivable [3] from the action of the classical string [74] 
(which has a linear potential and thus exerts a constant force), but follows more directly 
from symmetry arguments. 

Note that the junction motion need not be uniform. In the example above, one of 
the string pieces goes from the junction, via g to q'. At early times, the junction only 
experiences the pull of its immediate neighbour, g, and the direction of q' is irrelevant. 
However, as the gluon moves out from the origin, it loses energy to the string traced out 
behind it. From the point when all its energy has been converted to potential energy of the 
string and this information has propagated back to the junction, it will be the direction of 
q' which determines the direction of the 'pull' exerted by this string piece on the junction. 
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Figure 13: Example of colour assignments in a 7p collision with two interactions. Explicit 
colour labels are shown on each propagator line. In this example, the string system con- 
taining the junction is spanned by J connected to q^, to qvs, and via g to q', as can be seen 
by tracing each of the three colour lines to the junction. 



and not that of the gluon. In the general case, with arbitrarily many gluons, the junction 
will thus be 'jittering around', being pulled in different directions at different times. 

However, rather than trying to trace this jitter in detail — which at any rate is at or 
below what it is quantum mechanically meaningful to speak about — we choose to define an 
effective pull of each string on the junction, as if from a single parton with four-momentum 
given by [3]: 

n 

Ppull = y^,'Pi exp f-Z]^=l-^j/-^norm) , (66) 
i=l 

where the outermost sum runs over the parton chain which defines the string piece, from 
the junction outwards (in colour space), and where the sum inside the exponent runs over 
all gluons closer to the junction than the one considered (meaning it vanishes for i = 1). 
The energy normalization parameter -Enorm is by default associated with the characteristic 
energy stored in the string at the time of breaking, -Enorm — 1-5 GeV. Naturally, the energies 
Ej should be evaluated in the junction rest frame, yet since this is not known to begin with, 
we use an iterative sequence of successively improved guesses. 

With the motion of the junction determined, the fragmentation of the system as a 
whole may now be addressed. Since the string junction represents a localized topological 
feature of the gluon/string field, we would not expect the presence of the junction in the 
string topology to significantly affect the fragmentation in the regions close to the endpoint 
quarks. Specifically, in an event where each of the three endpoint quarks have large energies 
in the junction rest frame, the energies of the leading and hence hardest particles of each 
jet should agree, on average, with that of an equivalent jet in an ordinary two-jet event. 

The hadronization model developed in [3] ensures this by fragmenting each of the string 
pieces outwards-in, as for a normal qq string (in both cases opposite to the physical time 
ordering of the process). The leading quark of a string piece is combined with a newly 
created quark-antiquark pair to form a meson plus a new leftover quark, and so on. Parton 
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flavours and hadron spins are selected in a manner identical to that of the ordinary string, 
as are fragmentation functions and the handling of gluon kinks on the string pieces. 

However, junctions were not included in the original string model, so here a new pro- 
cedure needs to be introduced. If all three string pieces were fragmented in the above 
way until little energy was left in each to form more hadrons, then it would be extremely 
unlikely that the resulting leftover system of three unpaired quarks would just happen to 
have an invariant mass equal to that of any on-shell baryon. While one could in principle 
amend this by shuffling momentum and energy to other hadrons in the vicinity of the junc- 
tion, such a procedure would be arbitrary and result in an undesirable and large systematic 
distortion of the junction baryon spectrum. The way such systematic biases are avoided 
for ordinary qq strings is to alternate between fragmenting the system from the q end and 
from the q end in a random way, so that the hadron pair that is used to ensure overall 
energy-momentum conservation does not always sit at the same location. Thus, while the 
distortion is still local in each event, it is smeared out when considering a statistical sample 
of events. 

In the case of a junction system, such a procedure is not immediately applicable. Instead, 
we flrst fragment two of the three string pieces, from their respective endpoint quarks 
inwards. At the point where more energy has been used up for the fragmentation than is 
available in the piece, the last quark-ant iquark pair formed is rejected and the fragmentation 
is stopped. The two resulting unpaired quarks, one from each fragmented string piece, are 
then combined into a single diquark, which replaces the junction as the endpoint of the 
third string piece. Subsequently, this last string piece is fragmented in the normal way, 
with overall energy and momentum conservation ensured exactly as described for ordinary 
strings above. In order to minimize the systematics of the distortion and ensure that it is 
at all possible to produce at least two hadrons from this flnal string system, we choose to 
always select the highest energy string piece as the last to be fragmented. It was shown 
in [3, 77] that this asymmetry in the description does not lead to large systematic effects. 

In proton-proton collisions, two junction systems will be present, but it is physically 
impossible for these to be connected in colour. Hence, the hadronization of each system 
again proceeds exactly as described above. However, in pp collisions a new possibility 
arises, as depicted in Fig. 14. This simple example goes to illustrate that a junction and an 
'antijunction' may become colour-connected by the colour exchanges taking place in a given 
process. In such cases, the fragmentation of each of the junctions is no longer disconnected 
from what happens to the other one; instead the fragmentation of the system as a whole 
must be considered. The necessary generalization of the principles outlined above to the 
case of connected junction-junction systems [3] is not very complicated. 

As before, two of the three strings from a junction are fragmented flrst, outwards-in 
towards the junction, but in this case we always choose these two string pieces to be the 
ones not connected to the other junction. Diquarks are then formed around each junction 
exactly as before. What remains is a single string piece, spanned between a diquark at one 
end and an antidiquark at the other, which can be fragmented in the normal way. In fact, 
the only truly new question that arises at this point is how to generalize eq. (66) to describe 
the pull of one junction on another. Here gluons on the string between the two junctions 
are considered as normally, i.e. their momenta are added, with a suppression factor related 
to the energy of the intermediate gluons. The partons on the far side of the other junction 
also contribute their momenta, separately for each of the two strings, with an energy sum 
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Figure 14: Example of colour assignments in a pp collision, with explicit colour labels 
shown on each propagator line. Note that the blue colour line starting on the junction J is 
connected via the colour flow of the hard scattering to the antiblue colour line of J. 
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Figure 15: a) A string system (dashed lines) spanned between four quarks and containing 
a junction and an antijunction. b) The same parton configuration in colour space but with 
an alternative string topology. In a) the beam baryon numbers will still be present in the 
final state, while in b) they will have disappeared through annihilation. 



suppression now given by the intervening gluons on that particular string, plus the gluons 
on the junction-junction string. 

However, an alternative topology is also possible, where the junction and the antijunc- 
tion annihilate to produce two separate qq systems [3], as illustrated in Fig. 15b. While 
it is not clear from basic principles how often this should happen, it seems likely that, 
for a given event, the topology which has the minimal string length is the one selected 
dynamically. In this case, the string topology depicted in Fig. 15a would result when the 
qiq2 and q]^q2 opening angles are small, while the topology in Fig. 15b would result if the 
qiq^ and q2q2 opening angles are the small ones. Since, in the context we discuss here. 
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the quarks colour-connected to the junction will more often than not reside in the beam 
remnants, we do not expect annihilation between the incoming baryon numbers to be a 
large effect. Indeed, for the range of more realistic models that are investigated in Section 6 
below, junction-junction annihilation is a feature of less than 1% of the events at Tevatron 
energies. 

An ugly situation occurs in the rare events when the two junctions are connected by two 
colour lines. If these lines contain intermediate gluons, it would be possible but difficult to 
fragment the system, in particular when the energy of these gluons becomes small. Without 
intermediate gluons, a first guess would be that the junction and antijunction annihilate 
to give a simple string spanned between a quark and an antiquark endpoint, so that the 
original baryon numbers are lost. However, this assumes that the system starts out from a 
point in space and time, a commonly used approximation in the string language. Viewed in 
the transverse plane of the collision, the original positions of the junctions and of the hard 
scatterings involved could well be separated by distances up to a fm, i.e. the intervening 
strings could have energies up to a GeV. It may then be that the strings can break before the 
junctions annilate, so that a baryon-antibaryon pair nevertheless is produced. A detailed 
modeling would be required, beyond the scope of the current study, and possibly beyond 
the validity of the string framework, so for now we choose to reject these rare events. 

5.2 Initial-State Colour Correlations 

In the planar approximation, Nq ^ oo, a 2 ^ 2 process, such as gg gg, can receive 
contributions from several possible colour fiows, but the cross section for a colour fiow is 
uniquely defined, so that each fiow can be selected according to its appropriate weight [78]. 
Furthermore, in our leading-order parton showers, the colour fiow is well-defined in each 
branching. Within each separate interaction and its associated shower activity, the colour 
fiow can thus be selected unambiguously. In events with only one hard interaction, the 
colour of the shower initiator can also uniquely be hooked up to that of the beam remnant. 
Thus a knocked-out quark leaves behind a colour antitriplet beam remnant, a gluon leaves 
a colour octet beam remnant, and so forth. 

Unfortunately, once several interactions are allowed, there is no longer a unique answer 
how to hook up the different shower initiators with the beam-remnant partons. To illustrate, 
consider an incident meson out of which n gluons are kicked out. These gluons may be 
ordered in colour sequence such that the colour of the quark matches the anticolour of 
one gluon, which then has another colour that matches the anticolour of another gluon, 
and so on till the antiquark. Obviously there are n! such possible arrangement of the 
n gluons, each leading to a unique colour topology for the hadronizing partonic system. 
Perturbation theory has nothing to say about the relative probability for each of these 
configurations; the colour correlations we now consider arise at scales below the parton- 
shower cutoff Qq ^ 1 GeV. Further arrangements would exist if we allowed some of the 
above gluons to form a separate colour singlet, disconnected from the sequence between the 
quark and antiquark ends; some of these could contribute to diffractive topologies. 

So far we only considered the planar colour topologies of the Nc oo limit. The 
real-world Nq = 3 offers further complications. First, interference terms of order — 
modulated by kinematical factors — arise between different possible colour flows in hard 
processes. The more partons are involved, the more assignments need be made, and the 
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Figure 16: Example of how a given set of parton shower initiators could have been radiated 
off the initial baryon valence configuration, in the case of the 'purely random' correlations 
discussed in the text. In this example, the baryon number is disconnected from the beam 
remnant. Instead, it is the final-state partons connected to the colour lines of gi, g2, and 
g3 which determine how the junction moves and hence how the baryon number flows in the 
event. 



larger the total uncertainty. Even at the perturbative level it is thus no longer possible 
to speak of a unique colour arrangement. Second, the situation is more complicated for 
baryon beams. As described above, the initial state of a baryon, before any scatterings 
occur, is represented by three valence quarks, connected in a Y-shaped topology via a 
central junction which acts as a switchyard for the colour flow and which carries the net 
baryon number. This situation is illustrated in Fig. 12. Each of the gluons considered in 
the meson-beam example above may now be arranged in colour on either of the three string 
pieces, leading to a further multiplication of possibilities. 

We choose to address this question by determining a sequence of fictitious gluon emis- 
sions by which this configuration evolves (in colour space) to give rise to the parton shower 
initiators and beam remnant partons actually present in a given event. We here assume 
that only the minimal number of emissions required to obtain the given set of initiators and 
remnants is dynamically relevant. Further, since sea quarks together with their companion 
partners can pairwise be associated with a gluon branching below the parton shower cutoff, 
only gluon emissions remain to be considered. (This also means that a sea quark, in our 
model, can never form a colour singlet system together with its own companion.) 

5.2.1 Random Colour Correlations 

The simplest solution would be to assume that Nature arranges these correlations randomly, 
i.e. that gluons should be attached to the initial quark lines in a random order, see Fig. 16. 
In this case, the junction (and hence the baryon number) would rarely be colour connected 
directly to two valence quarks in the beam remnant, even in the quite common case that 
two such quarks are actually present (multiple valence quark interactions are rare). It 
should be clear that the migration of the baryon number depends sensitively upon which 
partons in the final state the junction ends up being connected to (see further Section 5.1). 
The conclusion is that if the connections are purely random, as above, then the baryon 
number will in general be disconnected from the beam remnant valence quarks. Hence 
the formation of a diquark in the beam remnant would be rare and the baryon number 
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Figure 17: p± spectra for junction baryons, a) with primordial k± switched on, and b) with 
primordial k± switched off. The shaded area represents the distribution in the new model 
of those junction baryons which arose by first forming a diquark in the beam remnant, 
cf. Section 4.2. 



of the initial state should quite often be able to migrate to small xp values, as previously 
illustrated in Fig. 11. One could expect this longitudinal migration to be accompanied by a 
migration in the transverse plane, such that the junction baryon should generally migrate 
to larger p± values when the junction is allowed to 'float' more. However, as Fig. 17a 
illustrates, no large differences in the total p± spectrum are apparent when comparing the 
new model (thick dashed) with the old Tune A (solid). 

The reason that no large transverse migration effect is visible, relative to the old model, 
is that the latter also has a broader leading-baryon p^-spectrum than for normal baryons, 
as follows. In the old model, the beam remnant diquark (around which the junction baryon 
forms in the fragmentation) always receives the full primordial kick from the hard 
interaction initiator, by default a Gaussian distribution with a width of 1 GeV. In the 
fragmentation process, when the baryon acquires a fraction of the diquark longitudinal 
momentum, it obtains the same transverse momentum fraction. Additional p± will be 
imparted to the junction baryon from the newly created quark, at the level of 0.36 GeV, 
but with some dependence on the momentum-space location of its nearest neighbour in 
colour space. Essentially, these two effects combine to yield the sohd curve in Fig. 17a. 

In the new model, we must distinguish between the old-model-like case when a diquark 
is formed in the beam remnant, on one hand, and those cases where the junction is 'free' 
to migrate, on the other. 

If a diquark is formed, then it consists of two undisturbed beam remnant quarks which 
are colour connected to the junction, and the situation is indeed very similar to the old 
model. The baryon forming around this diquark receives p± from three sources. Firstly, the 
diquark will have some intrinsic primordial k±, distributed according to eq. (45). Since the 
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diquark resides entirely within the beam remnant, this k± will always be at the level of the 
fragmentation pj_. Secondly, primordial k± will be imparted to the diquark by recoil effects 
from other beam remnant and initiator partons. In the case that primordial k± kicks are 
compensated for uniformly by all initiator and remnant partons, it is normally impossible for 
the diquark to acquire more than a fraction of the hardest interaction initiator's primordial 
k±. Even when k± compensation is more local, by straightforward combinatorics, the more 
initiators present in the event, the smaller the chances that the initiator parton(s) closest 
in colour to the diquark is associated with a scattering at large Q^. Hence, again according 
to eq. (45), it is apparent that the diquark usually will not receive a very hard primordial 
k± kick. Thus, such a diquark will in general have a smaller total primordial k± than a 
diquark in the old model. As before, the baryon will keep a large fraction of this diquark p± 
in the fragmentation process, as well as obtaining extra fragmentation p±. The net result 
is a softer junction baryon transverse momentum spectrum than in the old model, as can 
be verified by comparing the asymptotic slope of the shaded area in Fig. 17a with that 
of the solid curve. This conclusion is further established by the observation that, when 
primordial k± effects are not included, see Fig. 17b, indeed the spectrum of the old model 
becomes almost identical to that of the shaded region. In addition, it can already here be 
recognized that the junction baryon must have larger p_i_ in those events where a diquark is 
not formed, by comparing the slopes of the full junction baryon spectrum (dashed curves) 
with those of the shaded regions in either figure. We now study this further. 

If a diquark is not formed, then the junction may a priori be colour connected to partons 
going in widely different directions in the transverse plane. Nonetheless, as was described 
in Section 5.1, the fragmentation occurs in such a way that the junction baryon is always 
the last, i.e. normally slowest, hadron to be formed in either of the three directions. Hence, 
while the colour neighbours of the junction may themselves have large transverse momenta, 
this momentum will in general be taken by the leading hadrons formed in the fragmentation 
and not by the junction baryon. Unless two of those partons are going in roughly the same 
direction in if, the junction baryon itself will still only obtain a fairly small p±. The end 
result is a rather small p± enhancement, that is masked by the decreased primordial k±, 
Fig. 17. 

Thus, the main difference in the new model is that the beam baryon number can migrate 
longitudinally to a much larger extent than in the old model. Empirically, it may be 
desirable to be able to limit the degree to which this baryon number stopping occurs, 
and furthermore both perturbative and impact-parameter arguments allow much of the 
activity to be correlated in 'hot spot' regions that leave the rest of the proton largely 
unaffected. Therefore a free suppression parameter is introduced, such that further gluons 
more frequently connect to a string piece that has already been disturbed. In this way, 
gluons would preferentially be found on one of the three colour lines to the junction. This 
will reduce the amount of baryon number stopping and is an important first modification, 
but most likely it is not the only relevant ordering principle. 

5.2.2 Ordered Colour Correlations 

With the gluons connected preferentially along one of the three colour lines to the junction, 
we now address the question of their relative order along that line. If this order is random, 
then strings will in general be stretched criss-cross in the event. This is illustrated in 
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Figure 18: Example of initial-state colour correlations in an imagined event in which three 
gluons have been knocked out of an incoming hadron by colourless objects (for simplifica- 
tion), with no parton showers. In case a) the gluons have been randomly attached in colour 
to each other and to the beam remnant, as indicated by the dashed lines, whereas in case 
b ) the connections follow the rapidities of the hard scattering systems. 



Fig. 18a for a very simplified situation. However, it is unlikely that such a scenario catches 
all the relevant physics. More plausible is that, among all the possible final-state colour 
topologies, those that correspond to the smaller total string length are favoured, all other 
aspects being the same. One possible way of introducing such correlations is illustrated in 
Fig. 18b. 

In case a) of Fig. 18 there are four string pieces criss-crossing the rapidity range between 
the systems I2 and I3, while in case b) there are only two string pieces spanning this range. 
Hence the total string length should on average be smallest for the latter type of correlations. 
However, the rapidity distance is not the only variable determining the string lengths, also 
the transverse separations play a role. Moreover, when the interactions exchange colour 
between the colliding objects, there is no longer a unique correspondence between the colour 
flow of the initial state and that of the final state. 

To investigate these effects quantitatively, we consider three different possibilities for 
initial-state colour correlations (with the suppression of attachments breaking up the beam 
remnant applicable to all cases): 

1. Random correlations, as in Fig. 18a. 

2. Initiator gluons are attached preferentially in those places that order the hard scat- 
tering systems in rapidity, as in Fig. 18b. The rapidities are calculated at a stage 
before primordial k± is included. Hence, y = |ln^. 

For beam remnant partons, the rapidities are not yet known at the stage discussed 
here, since the initial-state colour connections are in our framework made before pri- 
mordial k± and beam remnant x values are assigned. However, beam remnant partons 
are almost by definition characterized by having large longitudinal and small trans- 
verse momenta. Thus, we assign a fixed, but otherwise arbitrary, large rapidity to 
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Figure 19: Example showing the colour flow produced by attaching the gluon g at the place 
indicated by the cross. 



each of the beam remnant partons, in the direction of its parent hadron. Finally, 
gluons are attached sequentially to the initial valence topology, with the attachments 
ordered by minimization of the measure 

= \yg-yi\ + \yg-y2\, (67) 

where Ug is the rapidity associated with the attached gluon and yi^2 are the rapidities 
associated with the partons it is inserted between. For those gluons which appear 
only as parents of sea quark pairs, the rapidity of the most central of the daughters 
is used. 

Note that, since the same hard-scattering rapidities are used for both beam remnants, 
the ordering in the two remnants will be closely correlated in this scenario, at least 
as long as only gluon-gluon interactions are considered. 

Initiator gluons are attached preferentially in those places that will give rise to the 
smaller string lengths in the flnal state. This is the most aggressive possibility, where 
the actual momentum separations of flnal-state partons, together with the full colour 
flow between the two sides of a hadronic collision, is used to determine which gluon 
attachment will result in the smallest increase in potential energy (string length) of 
the system, with each gluon being attached one after the other. The measure we use 
to deflne the increase in string length, for a particular attachment, is [3,81] 



AA = ln 



2 (Pci -Pci) (Pc2 - Pea) 



(68) 



where mo is a normalization constant, which drops out when comparing the string 
lengths of two different gluon attachments, and Ci (02) represents the flnal-state parton 
carrying the colour (anticolour) index of the attached gluon. To illustrate. Fig. 19 
shows the partons that enter the above expression for a speciflc example. Before 
the attachment, a single string piece is spanned between the final-state partons that 
carry the colour indices denoted ci and C2. After the attachment, there are two string 
pieces, one that is spanned between ci and ci, the other between C2 and C2, hence the 
increase in string length is given by the expression eq. (68). 
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As above, however, note that neither primordial k± nor beam remnant longitudinal 
momenta have yet been assigned at this stage. Simplified kinematics are therefore 
set up, to be used only for the purpose of determining the colour connections: the 
momentum remaining in the beam remnant on each side is divided evenly among the 
respective remnant partons (junctions are here treated simply as 'fictitious partons', 
receiving the same momentum as the 'real' remnant partons), and primordial k± 
effects are ignored. Thereby, parton pairs involving (at least) two partons in one of 
the beam remnants will come to have zero invariant mass, hence the total AA will 
be negative infinity for such pairs. Obviously, this is not desirable; one string piece 
with vanishing invariant mass should not affect the comparison, hence we impose 
a minimum invariant mass of rriQ for each string piece. If knowledge of the full 
kinematics of the final state was available a better choice could of course be made 
here. However, these two aspects are intertwined. The kinematics of the final state 
may depend on the colour connections assumed for the initial state (see Section 4.3 
above), and vice versa. Our choice has been to determine the initial-state colour 
connections first, and then subsequently construct the final-state kinematics, hence 
some approximation is necessary at this point. 

A variable which we have found to be sensitive to the colour connections in an event is 
the mean pj_ vs. charged multiplicity, (p_L)(?^ch) [!]• In scenarios with large string lengths, 
each additional interaction would result in a large increase in hadron multiplicity. This 
large multiplicity per interaction means that, in such scenarios, observed average charged 
multiplicities are reproduced with comparatively large values of p±o, i.e. with only a few 
parton-parton interactions taking place per event. Hence, correspondingly little perturba- 
tive p± is generated. On the other hand, in scenarios with smaller string lengths, compar- 
atively more interactions would be required to produce the same multiplicity, hence more 
perturbative p± would be generated per charged particle, bringing {p±){nch) up. 

In Fig. 20, we show the {p±) vs. rich distribution for each of the possibilities described 
above. In all cases, p±Q was first selected so as to give identical average multiplicities, 
corresponding to the multiplicity obtained with Tune A. Since the AA ordering results 
in the largest average number of interactions for a given multiplicity, it has the largest 
average p± per particle of the new scenarios, while the random ordering results in the 
smallest number of interactions. 

One also notes that Tune A, which more or less agrees with recent experimental data 
[56,58,61], shows an even steeper rise with rich than any of the new scenarios. This tune 
of the old model is such that the partons produced by subsequent scatterings will almost 
always be hooked up to the existing configuration in the way that minimizes the total string 
length. This is more or less like the AA ordering described above, but with the essential 
difference that the AA ordering only concerns the colour lines that are present in the initial 
state, while the ordering of parton attachments in the old model occurs in the final state, 
without any attempt at constructing a consistent colour flow in the event. 

From these observations, an interesting inference can be made. By the failure of even 
the AA ordering of the colour lines in the initial state to describe the {p±){nch) distribution, 
it appears that the colour flow in physical events cannot be correctly described by merely 
arranging the colour lines present in the initial state. We imagine two possible causes for 
this. Firstly, the initial-state showers associated with each scattering are constructed by 
backwards DGLAP evolution of each scattering initiator separately, down to the shower 
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Figure 20: {p±) vs. rich at the Tevatron for Tune A (solid lines), and for the new model 
with random (dashed lines), rapidity ordered (solid lines), and string length ordered (dotted 
lines) correlations in the initial state. Note that the origo of the plot is not at (0,0). For 
each of the new MI scenarios, p_Lo was selected to give the same average charged multiplicity 
as Tune A, with the same impact parameter dependence as Tune A (i.e. a double Gaussian 
matter distribution). 



cutoff scale. This does not take into account the possibility that the showers could be 
intertwined, i.e. that a parton at low virtuality, but above the shower cutoff scale, could 
have branched to give rise to two higher- virtuality scattering initiators. Secondly, one or 
more mechanisms causing colour exchanges between the showers may be active, both in the 
initial state as well as in the final state. Below we present some first studies related to the 
topic of colour exchanges, well aware that no simple solutions are to be expected. 

5.3 Final-State Colour (Re-) Connections 

To investigate how much more we need to 'mess around with the colours', we study a crude 
model of colour exchanges in the final state. Essentially, we rearrange the colour connections 
between the final-state partons in a manner that, taken to the extreme, will converge on 
that string configuration which has the smallest total 'string length', according to the A 
measure introduced above. This corresponds roughly to a minimization of the multiplicity 
produced when the parton system hadronizes. Note that we only apply this procedure to 
events where at least two interactions have occurred, to avoid the reconnections leading to 
large central rapdity gaps, i.e. diffractive topologies. 

For a given configuration of final-state partons in momentum and colour space, we apply 
an iterative procedure that successively brings down the total string length of the system, 
by the following steps: 
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1. First, we have assigned Les Houches Accord style colour tags [82] to all partons, so that 
each colour tag in the final state is matched by exactly one corresponding anticolour 
tag in the final state, with one string piece spanned between them. Junctions are 
special, since colour lines end there. In the following, we do not consider string pieces 
ending on junctions. 

2. Secondly, we decide on a fraction, F, of the colour tags present in the event for which 
we will attempt to make a reassignment. Note that F can be larger than one, since 
several different reassignments (n(n — 1)/2 for n colour tags, neglecting junctions) are 
normally possible. 

3. Next, we select two colour tags at random, ci and C2. Denoting the final-state parton 
carrying Ci colour (anticolour) by ii (ji) and the one carrying C2 colour (anticolour) 
by «2 (^2)5 we compute the combined string length. A, for the two string pieces ii — ji, 



By swapping e.g. the anticolours, a different string topology arises, ii — J2, ^2 — ji, 
with length 



If A' < A, the colour reassignment is accepted, otherwise the original assigments are 
kept. If the fraction of colour tags tried so far is smaller than F, a new pair of random 
colour tags is selected. 

4. Once the fraction F of colour tags has been tried, two things can happen. If at least 
one reconnection was made, then the colour topology now looks different, and the 
entire iteration is restarted. If no reconnection was made, the iteration ends. 

Briefly summarized, we thus introduce the fraction F as a free parameter that controls the 
strength of colour reconnections in the flnal state. 

As stated, this method is very crude and should not be interpreted as representing 
physics per se, but it does allow us to study whether a signiflcant effect can be achieved 
by manipulating the colour correlations to reduce the string lengths. As illustrated by 
the dashed histogram in Fig. 21, this is very much the case. Here, we have allowed a large 
amount of reconnections to occur, F = 1, adjusting p^o down so as to reproduce the average 
charged multiplicity of Tune A. As could be expected with this somewhat extreme choice of 
parameters, the new model now lies well above the data. (We shall return to more realistic 
tunings in Section 6.) 

However, this result should not be taken as evidence for the existence of colour reconnec- 
tions in physical events. Rather it allows us to infer that, by changing the colour structure 
of events, it should be possible to obtain agreement with the data within our framework. It 
is encouraging that, by studying and attempting to describe this distribution, we may learn 
interesting lessons concerning the highly non-trivial issue of colour flow in hadronic inter- 
actions. We therefore plan to go further, to construct more physically motivated models 
for colour rearrangement between partons, both in the initial state and in the flnal state, 
and also to allow for the possibility of intertwining the initial-state showers. 




(69) 




(70) 
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Tevatron: <p^>(N^J 

^ Old Ml: Tune A (p^^ = 2.0 GeV) 

V" New Ml: with FS reconnect (p^^ = 1 .90 GeV, F=1 .0) 

New Ml: no FS reconnect (p^g = 3.10 GeV, F=0.0) 

0.6 




Figure 21: {p±) vs. rich at the Tevatron for Tune A (solid line), and for the new model 
with (dashed line) and without (dotted line) final-state reconnections allowed. Both of the 
new models use rapidity ordering of the colour lines in the initial state and give the same 
average charged multiplicity as Tune A, with the same impact parameter dependence as 
Tune A. 



Although all of these issues appear almost hopelessly complicated from the point of 
view of pure QCD, the sahent features of the resulting physics may not be all that hard to 
penetrate. For instance, we imagine that the probability for two hard-scattering initiators 
to have originated from a common branching should be proportional to the probability 
that their spatial wavefunctions overlap, i.e. two very high-virtuality initial-state partons 
associated with different scatterings are most likely uncorrelated, since they only resolve 
very small distance scales in their parent hadron. On the other hand, two low-virtuality 
partons, even though associated with different scatterings, may very well have come from 
one and the same parent parton, since their wavefunctions are comparatively much larger. 
We will explore such and other ideas in a future study. 

6 Model Studies 

In this Section we concentrate on illustrating the properties of models that, as a baseline, 
roughly reproduce the charged multiplicity distribution of Tune A in pp collisions at a 
centre-of-mass energy of 1800 GeV. Our studies only concern inelastic nondiffractive events, 
i.e. essentially the same as the experimentally defined (trigger-dependent) "min-bias" event 
sample; we will here use the two concepts interchangeably. 

A general feature of the new multiple interactions modeling is that the added parton 
showers and the less efficient string energy minimization result in a higher multiplicity 
per interaction than Tune A. In order to arrive at the same average hadron multiplicity 
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Table 1: Parameters of the three models investigated in the text, and for Tune A where 
applicable. Also shown for each model is the average number of parton-parton interactions, 
('^int), the average number of final-state partons, (npart), and the average number of colour 
reconnections taking place, (^ircp), per min-bias collision at the Tevatron and at the LHC. 



as Tune A, without too much colour reconnection required in the final state, a generally 
larger p^o cutoff should be used. Table 1 lists three different tunes of the new framework, 
with successively smaller p±o values and with different schemes for the initial-state colour 
correlations. 

• The "Ran" model is based on a random ordering of the initial-state colour correlations, 
with a fairly large suppression of initiator gluon attachments to colour lines wholly 
within the beam remnant. Since only a minimal ordering of the colour correlations 
in the initial state are thus imposed, each additional interaction will ab initio give 
rise to a relatively large increase in hadron multiplicity. Therefore, a comparatively 
large cutoff p±o is used, and the F parameter — controlling the amount of final-state 
reconnections — is likewise chosen fairly large, so as to get the correct average charged 
multiplicity. 

• The "Rap" model uses the Ay measure introduced above to order the intial-state 
colour connections. p±o can thus here be slightly smaller, allowing more interactions 
on the average (with the same F fraction) for the same average charged multiplicity. 

• The "Lam" model employs the AA ordering of the initial-state colour correlations. In 
principle, this model should provide the most ordered initial state of the three, and 
thus allow a smaller p±Q and/or F. Unfortunately, the earlier-mentioned limitations, 
that the beam remnant kinematics is not fully fixed when the minimization is per- 
formed, leads to final string lengths which are not significantly shorter than for the 
Ay ordering. Choosing a smaller p±Q for this tune, the F fraction is consequently 
also required to be slightly higher, in order to reproduce the Tune A average charged 
multiplicity. 

Observe that all three models have a significant number of reconnections per event, cf. (nrcp) 
in Table 1. As a fraction of the total number of potential colour rearrangements it is below 
the 10% level, but one should keep in mind that several clusters of partons appear in 
reasonably coUimated jets, where reconnections would not be expected anyway. In this 
perspective, the amount of reconnections is quite significant. 

Common for the new models is a rather smooth overlap profile ExpOfPow(1.8), as 
compared to the more peaked double Gaussian of Tune A, this to better reproduce the 
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Figure 22: Multiplicity distributions for the Tevatron as obtained with Tune A (solid), and 
the Ran (dashed), Rap (dotted), and Lam (dash-dotted) models defined in Table 1. In all 
cases, the average charged multiplicity is 49.5 (within ±0.5). 



shape of the Tune A multiplicity distribution. In addition, all of the new models assume 
that primordial k± kicks are compensated uniformly among all other initiator partons, that 
composite objects can only be formed in the beam remnant by valence quarks, and that 
initial-state colour connections breaking up the beam remnant are suppressed, so that the 
relative probability of attaching a gluon between two remnant partons (i.e. breaking up the 
remnant) as compared to an attachment where at least one 'leg' is outside the remnant is 
0.01, whenever the latter type of attachment is possible. 

Fig. 22 shows the Tevatron multiplicity distributions of these models, as compared to 
Tune A. It is apparent that, while the average charged multiplicity is the same, the shape 
of the Tune A multiplicity distribution is not exactly reproduced by any of the models here 
investigated. This should not be taken too seriously; our aim is not to present fuU-fiedged 
tunes, rather it is to explore the general properties of the new framework, and how these 
compare with those of Tune A. 

Below, we first present comparisons for pp min-bias collisions at 1.8 TeV CM energy, 
highlighting the differences (and similarities) between the new models and Tune A. There- 
after, we apply the same models to the case of pp min-bias events at 14 TeV CM energy. 
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Figure 23: The number of jets as a function of jet Ej_ (for Ej_ > 5GeV) in min-bias collisions 
at the Tevatron. Results are shown for Tune A (solid), Ran (dashed), Rap (dotted), and 
Lam (dash-dotted), as defined in Table 1. 



6.1 Comparisons at the Tevatron 

Despite the differences between the old and new frameworks, the starting point in both 
cases is still that of a perturbative sequence of p^-ordered scatterings. Especially for the 
hardest partons there should thus be next to no difference between the old and new frame- 
works. An illustration of this is given in Fig. 23, where the probability of finding a jet 
with transverse energy Ej_ is plotted against E_i_. A simple cone algorithm with cone 
size AR = ^At]^ + A</)2 = 0.7 has been used to cluster the jets, and only particles with 
|?7| < 2.5 are included. As can be observed, there is hardly any difference between Tune 
A and the new models here. Further, the models exhibit similar charged particle spectra 
both in transverse momentum and in rapidity. Fig. 24, with a slightly harder p± spectrum 
and, consequently, a slightly more central y one in the new models. 

Having thus convinced ourselves that the overall features of the models agree, we turn 
to the aspects in which the new and old scenarios are expected to differ. One significant 
change is the possibility to knock out several valence quarks from the beam hadron. To 
quantify, the number of quarks (excluding diquarks) in the final state is illustrated in Fig. 
25a. As can readily be observed, final-state quarks are much more abundant in the new 
scenarios. Valence quark interactions are, however, not the only cause of this. Diquarks are 
not included in Fig. 25a, and (as discussed in Section 4.2) these are less frequently formed in 
the new scenarios, hence more of the quark content here appears as individual quarks in the 
final state. Further illustration of this is given by Fig. 25b, where the x = IE j ^ values of 
all final-state partons, including diquarks, are shown. The peak towards low values comes 
mainly from gluons and is the same for Tune A and the new models. However, some of the 
content of the x = 1 peak in Tune A has vanished in the new models and has been replaced 
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a) b) 

Figure 24: Charged particle p± and y spectra at the Tevatron. Results are shown for Tune 
A (solid), Ran (dashed), Rap (dotted), and Lam (dash-dotted), as defined in Table 1. 
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Figure 25: a) Number of final-state quarks (not counting diquarks) and b) final-state parton 
X values at the Tevatron. Results are shown for Tune A (solid). Ran (dashed). Rap (dotted), 
and Lam (dash-dotted), as defined in Table 1. 
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Tevatron: d(p+n-pbar-nbar)/dy Tevatron: y - Junction Baryons 




a) b) 

Figure 26: a) n{p) + n(n) — n(p) — n(n) and b) junction baryon rapidity distributions at 
the Tevatron. Results are shown for Tune A (sohd), Ran (dashed), Rap (dotted), and Lam 
(dash-dotted), as defined in Table 1. 



by a larger plateau at intermediate x, since a fraction of large-x diquarks has been split up 
into individual valence-like quarks. 

Another point where the models differ is in the treatment of the beam remnant, espe- 
cially concerning the flow of baryon number. Fig. 26a shows the distribution of baryons 
minus the distribution of antibaryons as a function of rapidity, illustrating one way of exper- 
imentally probing the location of the beam baryon numbers in the flnal state. As expected, 
the distribution is more peaked in the old scenario, where the beam baryon number is 
'locked' inside the escaping remnant diquark. Further illustration of the baryon number 
stopping is given by Fig. 26b, which shows the rapidity distribution of the flnal-state baryon 
carrying the baryon number of the incoming proton. Note especially the long tail in the 
new models, even extending to negative rapidities. The height of this tail depends on a set 
of model parameters, such as the choice of initial- and flnal-state colour correlations and 
the rules for diquark formation, cf. Fig. 11 fore more extreme scenarios. 

The fact that the baryon number is not so closely associated with the valence quarks in 
the new model has an interesting side effect: apart from migrating in physical space, the 
initial-state baryon number can also 'migrate' in flavour to a much larger extent than before. 
When the junction is resolved, the flavour selection rules of ordinary string fragmentation 
take over in determining the flavour composition of the junction baryon. Thereby, it be- 
comes possible for the junction baryon to have a larger strangeness number S than before, 
as shown in Fig. 27. The last bin, 15*1 = 3, is actually empty for Tune A, since in the old 
model it is impossible to produce an Q~ from the incoming beam baryon number. We note 
that heavy- ion experiments do observe a ratio > 1 [83], which is in contradiction 

with the old string model but which would be more in line with expectations based on the 
junction scenario introduced here. 
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Figure 27: Junction baryon strangeness. Results are shown for Tune A (solid), Ran 
(dashed), Rap (dotted), and Lam (dash-dotted), as defined in Table 1. 
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Figure 28: Average p_[_ as a function of charged multiplicity for min-bias collisions at the 
Tevatron. Results are shown for Tune A (solid). Ran (dashed). Rap (dotted), and Lam 
(dash-dotted), as defined in Table 1. 
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Finally, Fig. 28 shows the (p±)(?T-ch) distributions. As previously noted, Fig. 21, this 
distribution is very sensitive to the colour correlations present at the hadronization stage. 
Since this is one of the major open issues remaining, it is not surprising that the agreement 
here is far from perfect: the new models exhibit a too early rise to a too low plateau, as 
compared to Tune A. While it is possible to obtain a much better agreement by varying the 
scheme adopted for the colour reconnections in the final state, our attempts in this direc- 
tion have so far led to poorer descriptions of other distributions, the charged multiplicity 
distribution in particular. Moreover, the colour reconnection scheme adopted here is meant 
only as an instructive example, not as a model of the physics taking place. Our plan is to 
continue the study of colour correlations in more depth. In the context of such studies, it is 
encouraging to note that the acute sensitivity of the {p±){nch) distribution to these aspects 
makes its proper description a prime test for any physical model of the colour flow. 

6.2 Comparisons at the LHC 

Turning now to the situation at the LHC, the lack of experimental constraints increases the 
uncertainties. We here focus on only a subset of these, assuming the same energy scaling of 
the p±Q cutoff as for Tune A, proportional to -Ecm^; using the same parton distributions 
in all cases. We note that these aspects do constitute important sources of uncertainty in 
our ability to make trustworthy 'forecasts' for the LHC. Here, however, our aim is merely 
to compare alternative scenarios under identical boundary conditions. 

Fig. 29 shows the LHC charged multiplicity distribution for the same models as used 
for the Tevatron, but with the p±Q cutoff scaled to the LHC CM energy. The new models 
exhibit average multiplicities of 6-8 more charged particles per event than the Tune A 
value, {rich) = 81. This illustrates a general effect in the new models, which is due to the 
increased shower activity arising from associating also the sub-leading interactions with 
initial- and final-state cascades. The larger available phase space at higher energies implies 
that showers are more important at LHC, cf. the average number of final-state partons, 
(ripart), in Table 1. All else being equal this causes the multiplicity to increase more rapidly 
with energy than in the old model. (The strange bump on the Tune A distribution at low 
multiplicities is merely an artifact of the way parton distributions at low are handled in 
that model.) 

The addition of parton showers also increases the total amount of partonic transverse 
energy, but owing to a partial cancellation of the effects of radiation in the initial state 
(boosting some partons to larger p±) and in the final state (jet broadening), the jet rates 
come out similar, as depicted in Fig. 30 (adopting the same cone algorithm and \rj\ < 2.5 
region as before). 

Fig. 31 compares the junction baryon pj_ distributions at the Tevatron (left plot) and at 
the LHC (right plot), for Tune A and the new models. An interesting difference is that the 
junction baryon can be significantly harder in p_i_ at the LHC than at the Tevatron in the 
new models, whereas Tune A exhibits spectra which are almost identical between the two 
energies. This is due to the intrinsic difference between the way primordial k± is treated 
in the two frameworks. In the old model, the width of the primordial k± distribution for 
the parton initiating the hardest scattering is fixed, to 1 GeV by default, hence there is no 
mechanism that would allow the junction baryon spectrum to depend on the CM energy (at 
sufficiently high energies that energy-momentum conservation effects can be neglected) . In 
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Figure 29: Multiplicity distributions for the LHC as obtained with Tune A (solid), and the 
Ran (dashed), Rap (dotted), and Lam (dash-dotted) models defined in Table 1. For Tune 
A, the average charged multiplicity is {rich) = 81, whereas for the new models it is in the 
range 87-89. 



the new model, the amount of primordial k± given to initiators depends on the of their 
associated hard scattering. With the increased phase space at the LHC, more primordial 
kj_ is thus imparted by recoil effects to the junction baryon than at the Tevatron, hence the 
p± spectrum becomes harder. 

Also the junction baryon longitudinal migration shows some difference. Comparing the 
Tevatron junction baryon rapidity distribution. Fig. 26b above, with the LHC one. Fig. 32, 
we may distinguish two components. One is the peak at large rapidities, which corresponds 
to an (effective) diquark fragmentation and which is only shifted outwards in rapidity 
relative to the Tevatron by the increased energy. The other is the tail to central rapdities, 
which corresponds to baryon stopping. This tail does increase with energy, following the 
increase in the average number of interactions. 

Finally, we show the {p±){nch) distributions in Fig. 33. The same qualitative behaviour 
as at the Tevatron is apparent: the new models exhibit an earlier rise to a lower plateau, 
as compared to Tune A. Again, it is premature to draw any strong conclusions, in view of 
the still simple-minded description of the colour flow that we have included here. Further 
and more detailed studies of possible colour correlation mechanisms in hadronic collisions 
will be required in order to fully understand these aspects. 
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Figure 30: The number of jets as a function of jet E± (for E± > 5GeV) in min-bias collisions 
at the LHC Results are shown for Tune A (solid), Ran (dashed), Rap (dotted), and Lam 
(dash-dotted), as defined in Table 1. 
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Figure 31: Junction baryon -p^ spectrum at the Tevatron and h) the LHC. Results are 
shown for Tune A (solid). Ran (dashed). Rap (dotted), and Lam (dash-dotted), as defined 
in Table 1. 
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Figure 32: Junction baryon rapidity distributions at the LHC. Note: at the LHC both 
beam baryon numbers are included in the figure, whereas in the Tevatron plots, Fig. 26, 
the antibaryon number is not. Results are shown for Tune A (solid). Ran (dashed). Rap 
(dotted), and Lam (dash-dotted), as defined in Table 1. 
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Figure 33: Average p±_ as a function of charged multiplicity for min-bias collisions at the 
LHC. Results are shown for Tune A (solid). Ran (dashed). Rap (dotted), and Lam (dash- 
dotted), as defined in Table 1. 
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7 Conclusion and Outlook 



Only in the last few years have multiple interactions gone from being a scientific curiosity, 
by most assumed relevant only for some rare topologies of four-jet events, to being ac- 
cepted as the key element for understanding the structure of underlying events. However, 
this leaves a lot of questions to be addressed, such as: 

(i) What is the detailed mechanism and functional form of the dampening of the pertur- 
bative cross section at small p_L? 

(a) What is the energy dependence of the mechanism(s) involved? 

(in) How is the internal structure of the proton reflected in an impact-parameter-dependent 
multiple interactions rate, as manifested e.g. in jet pedestal effects? 

(iv) How can the set of colliding partons from a hadron be described in terms of correlated 
multiparton distribution functions of flavours and longitudinal momenta? 

(v) How does a set of initial partons at some low perturbative cutoff scale evolve into such 
a set of colliding partons? Is standard DGLAP evolution sufficient, or must BFKL/CCFM 
effects be taken into account? 

(vi) How would the set of initiators correlate with the flavour content of, and the longitu- 
dinal momentum sharing inside, the left-behind beam remnant? 

(vii) How are the initiator and remnant partons correlated by confinement effects, e.g. in 
primordial k±l 

(via) How are all produced partons, both the interacting and the beam- remnant ones, cor- 
related in colour? Is the large number-of-colours limit relevant, wherein partons can be 
hooked up into strings representing a linear confinement force? 

(ix) How is the original baryon number of an incoming proton reflected in the colour topol- 
ogy? 

(x) To what extent would a framework with independently fragmenting string systems, as 
defined from the colour topology, be modified by the space-time overlap of several strings? 

Tentative answers to some of the questions are provided by the Tune A of the Pythia 
multiple interactions framework. Thus we now believe that: 

• The matter overlap when two hadrons collide can be described by an impact-parameter 
dependence more spiked than a Gaussian but less so than an exponential. 

• The p_Lo regularization scale does increase with energy. 

• The colours of final-state partons are not random but correlated, somehow, to give a 
reduced string length. 

This still leaves many questions unanswered. Worse, existing event generators would not 
even address many of the relevant issues, at least not in a deliberate or realistic fashion. In 
this article we have therefore tried to take the next step towards a better understanding of 
the structure of a hadronic event, addressing several of the points above. This in particular 
has concerned the correlations between initiator and remnant partons in the hadron beams, 
in terms of flavour, longitudinal and transverse momenta. Colour correlations have also 
been studied, and here it appears that the final-state partons need be involved as well. The 
complexity of the colour issues is tremendous, however, and we do not consider the studies 
finished in this area. 

A specific new topic addressed is that of baryon number flow. Data from hadronic col- 
lisions, and even more from heavy-ion ones, show large excesses of baryon over antibaryon 
production in the central rapidity region of events [84], suggesting a significant influx of 
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baryon number from the high-rapidity coUiding beams, more than would be expected from 
standard quark/diquark fragmentation models. When the junction is introduced as a topo- 
logical feature of the colour field in the baryon, however, the fate of the baryon number 
of an incoming beam particle may partly or wholly decouple from that of the valence 
quarks [1,85]. We have here demonstrated that the junction topology in combination with 
multiple interactions can induce quite large rapidity shifts, of the desired kind. 

The problem may actually be the opposite, i.e. not to move the baryon number by 
too much. To this end, we have assumed a suppression of interactions that affect several 
of the three colour chains that connect the valence quarks to the junction. This could 
be an impact-parameter-related effect, that not the whole proton is involved in the hard 
processes. If so, the suppression should be less pronounced in heavy-ion collisions, where the 
interactions of a proton with several nucleons in the other nucleus could occur at different 
positions in the transverse plane and thereby affect different chains. Obviously it would be 
a major undertaking to construct a complete model for heavy-ion collisions to study these 
ideas, but we hope in the future to be able to present a simple study of the baryon number 
flow. 

Another open issue is that of intertwined initial-state showers, whereby two seemingly 
unrelated partons, each undergoing a hard scattering, reconstruct back to come from a 
common shower ancestor. With the new p^-ordered showers now being implemented in 
Pythia [86] we intend to introduce enough flexibility that such issues could be addressed. 

This will also further constrain the initial-state colour flow. The possibility of final-state 
colour reconnections remains, however, and has been proposed as a mechanism to introduce 
diffractive topologies in a number of processes [87]. One here needs to better understand 
how much reconnections are allowed/required, and of what character. 

We see that much work remains, before the physics of the underlying event is truly 
understood. Progress will not be possible without a constructive dialogue between theory 
and experiment. We have frequently had reason to mention Tune A as a role model here, 
because it offers a convenient reference that more sophisticated models can be tested against, 
without the need to know the details of the CDF detector. However, only a few distributions 
went into the tune, and so we do not know what to aim for in many other respects. 

To give one specific example, it would be valuable to have information on the 'lumpiness' 
of the underlying event, such as n-jet rates as a function of some jet resolution parameter, 
similarly to e"'"e~-annihilation QCD analyses. One would there hope for an intermediate 
resolution region, between the coarse one that is dominated by the perturbative QCD 
structure and the fine one that mainly is sensitive to hadronization details, where the 
structure of the multiple interactions would play a key role. An understanding of this 
lumpiness is related to the fluctuations in the jet pedestal, and thereby to the smearing of 
jet energies in SUSY searches, say. It all hangs together . . . 

In summary, striving for a better understanding of the physics of the underlying event 
is both interesting and useful. Interesting because it forces us to consider many issues 
normally swept under the carpet, and to confront dramatically different scenarios. Useful 
because it ties in with so many other physics analyses at hadron colliders. So there is plenty 
of interesting and useful work ahead of us before the picture has clarified completely! 
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